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R.  Saeks  and  R.-w,  Liu 
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Fault  Diagnosis  in  Electronic  Circuits 


R.  Saeks  R.-w.  Liu 

Department  of  Electrical  Engineering  Department  of  Electrical  Engineering 
Texas  Tech  University  Notre  Dame  University 

Lubbock,  Texas  79409  Notre  Dame,  IN.  46556 

During  the  past  quarter  century  the  engineering  community  has  been 
witness  to  tremendous  strides  in  the  art  of  electronics  design.  The  graph¬ 
ical  algorithms  of  the  previous  generation  have  given  way  to  the  modern 
CAD  package,  the  breadboard  has  been  subsumed  by  the  simulator.  Indeed, 
even  the  universal  building  block  has  become  a  reality.  To  the  contrary 
electronics  maintenance  has  changed  little  since  the  day  of  the  vacumm 
tube,  remaining  the  responsibility  of  the  experienced  technician  with 
scope  and  multimeter.  As  such,  our  ability  to  design  a  complex  electronic 
circuit  is  quickly  out-distancing  our  ability  to  maintain  it.  In  turn, 
the  price  reductions  which  have  accompanied  modern  electronics  technology 
have  been  paralleled  by  increasing  maintenance  and  operations  costs.  In¬ 
deed,  many  industries  are  finding  that  the  life  cycle  maintenance  costs 
for  their  electronic  equipment  now  exceed  their  original  capitol  invest¬ 
ment. 

Given  the  above,  it  is  quickly  becoming  apparent  that  the  electronics 
maintenance  process,  like  the  design  process,  must  be  automated.  Un¬ 
fortunately,  the  50  years  of  progress  in  circuit  theory,  on  which  our 
electronics  design  automation  has  been  predicated,  does  not  exist  in  the 
maintenance  area.  As  such,  the  past  decade  has  witnessed  the  inauouration 
of  a  basic  research  program  to  lay  the  foundations  for  a  theory  of  electronic 
maintenance  and  a  parallel  effort  to  develop  operational  electronic 
maintenance  codes. 
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Thus  far  the  greatest  success  has  been  achieved  in  the  digital  electronic 

area,  wherein  the  finite  state  nature  of  the  UUT  (unit  under  test)  may  be 

exploited.  Typically,  one  assumes  that  all  failures  manifest  themselves  in 

the  form  of  component  outputs  which  are  either  "stuck-at-one"  or  "stuck-at- 

c 

zero"  and/or  shorts  and  opens.  Under  such  an  assumption  a  theory  for 
digital  system  maintenance  has  been  developed  and  practical  fault  diagnosis 
algorithms  are  in  the  formative  stages  of  development.  Typically,  one 
hypothesizes  some  limit  on  the  number  of  simultaneous  faults  and  then  simu¬ 
lates  the  responses  of  the  UUT  to  a  family  of  test  vectors  for  each  allowed 
combination  of  faults.  The  actual  responses  of  the  UUT  are  then  compared  with 
the  simulated  responses  to  locate  the  failure.  Although  lacking  in  asthetic 
appeal  the  above  approach,  termed  fault  simulation,  is  ideally  suited  for 
the  maintenance  environment,  wherein,  the  actual  simulation  process  need  only 
be  done  once  at  the  factory  or  a  maintenance  depot  with  the  simulated 
response  data  being  distributed  via  magnetic  tape  to  the  various  field  loca¬ 
tions  where  the  actual  test  is  conducted.  As  such,  with  the  aid  of  some 
sophisticated  software  engineering,  this  apparently  "brute  force"  approach 

4 

to  the  fault  diagnosis  problem  has  slowly  evolved  into  a  workable  concept. 
Indeed,  at  the  present  time  a  number  of  automatic  test  program  generators 
which  classify  faults,  choose  test  vectors,  and  carry  out  the  appropriate 
simulation  (often  in  a  parallel  processing  mode),  are  commercially  available 
and,  as  such,  the  automated  maintenance  of  digital  electronic  circuits  is 
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becoming  a  reality. 

Unfortunately,  the  above  described  success  in  the  digital  world  has 
not  been  paralleled  by  progress  in  the  analog  world.  Indeed,  test  engineers 
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complain  that  while  80%  of  the  boards  are  digital,  80%  of  their  headaches 
are  analog  and  hybrid.  This  difficulty  arises  from  a  number  of  character¬ 
istics  of  the  analog  problem  which  are  not  encountered  in  digital  circuits. 

«  Indeed,  in  an  analog  circuit: 

(i)  there  is  a  continuum  of  possible  failures, 

(ii)  a  component  may  be  "in  tolerance"  but  not  nominal, 

»  (iii)  complex  feedback  structures  are  encountered, 

(iv)  simulation  is  slow  and  costly, 

(v)  post-fault  component  characteristics  may  not  be  known, 

•  (vi)  and  a  fault  in  one  component  may  induce  an  apparent  fault  in 

another. 

Items  (i)  and  (ii)  imply  that  an  extremely  large  number  of  simulations  will 

*  be  required  for  analog  testing.  Items  (iii)  and  (iv)  suggest  that  these 
simulations  will  be  far  more  expensive  than  similar  digital  simulations. 
Finally,  items  (v)  and  (vi)  indicate  that  the  simulation  of  a  post-fault 
circuit  by  itself  may  not  be  a  tractable  problem.  As  such,  it  is  by  no 
means  clear  that  the  kind  of  "brute  force"  fault  simulation  algorithm 
associated  with  the  digital  problem  will  be  applicable  to  the  analog  or 
hybrid  case. 

As  an  alternative  to  fault  simulation,  a  number  of  academic  researchers 

have  proposed  a  variety  of  "post  test"  fault  diagnosis  algorithms,  wherein, 

an  "equation  solving  like"  algorithm  is  used  to  locate  the  faulty  component 

2  8 

given  the  test  data  from  UUT.  ’  Although  these  algorithms  are,  in  some 
sense,  "smarter"  than  the  simulation  algorithms,  most  of  the  required  com¬ 
puting  must  be  done  in  the  field  after  the  UUT  has  been  tested.  Moreover, 
these  computational  requirements  must  be  replicated  each  time  a  unit  fails. 
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As  such,  the  success  of  such  "post  test"  algorithms  is  contingent  on  re¬ 
ducing  their  computational  requirements  to  a  bare  minimum.  Although  no 
system  is  yet  operational,  with  the  aid  of  the  powerful  linear  circuit 
theory  developed  over  the  past  half  century,  a  computationally  efficient 

solution  to  the  fault  diagnosis  problem  for  linear  analog  circuits  appears 
1  2 

to  be  within  reach.  ’  Unfortunately,  no  such  light  exists  at  the  end  of 
the  nonlinear  tunnel,  wherein  progress  appears  to  be  limited  by  a 
"computational  complexity/test  point"  bound. 

Not  surprisingly,  the  computational  cost  of  an  analog  fault  diagnosis 

algorithm  is  an  inverse  function  of  the  number  of  test  points  at  which 

measurements  of  the  UUT  may  be  made.  Indeed,  if  one  lets  n  be  a  measure  of 

UUT  complexity  (which  may  loosely  be  taken  to  be  the  total  number  of 

terminals  for  all  of  the  circuit  components),  then  if  one  has  access  to 

0(n)(t1)  test  points  the  fault  diagnosis  problem  can  be  resolved  using 

linear  algorithms. Moreover,  by  combining  such  algorithms  with  the 

above  mentioned  linear  algorithms,  acceptable  computational  efficiency  can 

be  obtained  with  0(m)  test  points  where  m  is  a  measure  of  the  complexity 

of  the  "nonlinear  subsystem"  of  the  UUT.^’^  Although  such  algorithms  can 

be  effective  on  the  typical  academic  example  a  "real  world"  PC  (printed 

circuit)  board  does  not  iave  terminal  space  for  the  20  or  30  test  points 

which  are  required  even  for  a  routine  board  made  up  of  discrete  components 

and/or  SSI  (Small  Scale  Integration).  Although  the  problem  can  be  partially 

alleviated  by  making  internal  measurements  with  the  aid  of  a  "bed-of-nails” 

tester  it  has  been  our  experience  that  such  testers  cause  as  many  failures 

as  they  locate  while  their  applicability  to  two-sided,  multilayer,  and 

t  f(n)  =  0(n)  means  f  increases  in  the  order  of  n;  more  precisely, 

|f(n)|  <_ c { n j  for  some  c  >  0. 
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coated  boards  is  severly  limited.  As  such,  we  would  like  to  limit  the 
number  of  test  points  required  by  an  analog  fault  diagnosis  algorithm 
should  inc  ease  at  a  rate  of  no  greater  than  0(n^2).  A  further  study  of 
the  possible  tradeoff  between  test  points  and  computational  cost  appears 
in  references  11  and  12. 

Unfortunately,  all  computationally  acceptable  "post  test"  algorithms 
which  have  thus  far  been  proposed  have  test  point  requirements  which  grow 
linearly  with  UUT  complexity  (assuming  that  m  grows  linearly  with  n).  As 
such,  many  researchers  are  looking  at  the  classical  fault  simulation  al¬ 
gorithms  with  renewed  vigor.  Indeed,  these  algorithms  have  minimal  on¬ 
line  computational  costs,  while  the  number  of  test  points  employed  can 

2/1 

easily  be  kept  below  0(n  ).  The  difficulty  lies  with  the  required  number 

of  simulations  and  the  development  of  decision  algorithms  which  will  allow 
us  to  "interpolate"  between  simulated  data  points. 

Thus,  while  the  state-of-the-art  in  digital  diagnosis  is  fast  maturing, 
a  serious  investigation  of  analog  fault  diagnosis  problems  is  only  just 
beginning.  Indeed,  a  satisfactory  fault  diagnosis  code  for  linear  analog 
circuits  has  yet  to  be  demonstrated  while  the  nonlinear  problem  has  yet  to 
progress  beyond  the  basic  research  stage. 
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FAULT  DIAGNOSIS  -  A  NONLINEAR  SYSTEMS  PROBLEM 


R.  Saeks 


Fault  Diagnosis  -  A  Nonlinear  Systems  Problem 


R.  Saeks 

Department  of  Electrical  Engineering 
Texas  Tech  University 
Lubbock,  Texas  79409 


Summary 

Conceptually,  the  fault  analysis  problem  for  an  analog  circuit  or 
system  amounts  to  the  measurement  of  a  set  of  externally  accessible  para¬ 
meters  of  the  system  from  which  one  desires  to  determine  the  internal  system 

★ 

parameters  or  equivalently  locate  the  failed  components  as  illustrated  in 
Figure  1 .  Here,  the 


Figure  1.  Conceptual  Model  of  Fault  Diagnosis  Problem. 

measurements,  m^ ,  may  represent  data  taken  at  distinct  test  points  or  alter¬ 
natively,  data  taken  at  a  fixed  test  point  under  different  stimuli. 

Since  the  problem  of  determining  the  values  of  the  failed  components  is  usually 
straightforward,  once  the  failures  have  been  located,  the  exact  determination  of 
all  internal  component  parameters  is  essentially  equivalent  to  the  problem  of 
"simply”  locating  the  failed  components. 
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Similarly,  the  r..  represent  parameters  characterizing  the  various  internal 
system  components.  Here,  a  single  parameter  may  characterize  an  entire 
component,  say  a  resistance,  capacitance  or  inductance.  Alternatively,  a 
component  may  be  represented  by  several  parameters:  the  h-parameters  of  a 
transistor,  the  poles  and  gain  of  an  op-amp,  etc.  In  general,  one  models 
a  system  component  by  the  minimum  number  of  parameters  which  will  allow 
the  failure  to  be  isolated  up  to  a  "shop  replaceable  assembly"  with  all 
"allowed"  system  failures  manifesting  themselves  in  the  form  of  some  para¬ 
meter  change. 

To  solve  the  fault  diagnosis  problem,  one  then  measures  m  =  col(m^) 
and  solves  a  nonlinear  algebraic  equation 

1 .  m  -  F(r) 

for  r  =  col(r^)  to  diagnose  the  fault.  Note,  the  function,  F,  is  nonlinear 
even  for  linear  systems,  however,  for  linear  time-invariant  systems  the 
function,  F,  can  be  expressed  analytically.  More  generally,  in  the  nonlinear 
case,  one  can  evaluate  F{r)  for  any  given  parameter  vector,  r,  with  a  simu¬ 
lator,  and  thus  solve  1.  numerically,  even  though  F  has  no  analytic  expression. 

Although  one  does  not  usually  formulate  the  fault  diagnosis  problem  in 
terms  of  the  above  described  equation  solving  notation,  this  formulation  is 
equivalent  to  the  classical  fault  simulation  concept.  Indeed,  fault  simula¬ 
tion  is  simply  a  search  algorithm  for  solving  1.  Here,  one  precomputes 
m  =  f(r)  for  each  allowable  faulty  parameter  vector  r  and  then  compares  the 
measured  m  with  the  simulated  m's,  stored  in  a  fault  dictionary,  to  solve 
equation  1 . 

#  * 

By  allowable  faults  we  mean  all  possible  parameter  vectors,  r,  which  satisfy  a 

specified  set  of  fault  hypotheses.  These  typically  restrict  the  maximum  number 

of  component  parameters  which  are  simultaneously  out  of  tolerance  and  the  type 

of  failure  (open  circuit,  short  circuit,  small  change,  etc.) 
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Although  the  above  described  approach  to  fault  simulation  has  been 
★  ★ 

successful  when  applied  to  digital  system,  there  is  considerable  question 
surrounding  its  applicability  to  analog  circuits  and  systems.  The  problem 
is  two-fold.  First,  rather  than  simply  failing  as  a  one  or  zero,  an  analog 
parameter  has  a  continuum  of  possible  failures.  Secondly,  unlike  a  digital 
system  wherein  a  component  is  either  good  or  bad,  in  an  analog  system,  a  com¬ 
ponent  parameter  is  either  in  tolerance  or  out  of  tolerance.  As  such,  for 
each  hypothesized  failure,  it  may  prove  necessary  to  do  an  entire  family 
of  Monte  Carlo  simulations  in  which  the  values  of  the  good  components  are 
randomly  chosen  within  their  tolerance  limits.  Although,  at  the  present 
time  we  have  insufficient  practical  experience  to  determine  the  precise 
number  of  fault  simulations  required  for  analog  fault  diagnosis,  it  is  esti¬ 
mated  that  the  number  of  simulations  required  for  an  analog  system  will  ex¬ 
ceed  the  number  of  simulations  required  for  a  digital  system  of  similar 
complexity  by  a  factor  ranging  between  two  and  six  order  of  magnitude.  As 
such,  the  fault  simulation  concept  which  has  proven  to  be  so  successful  for 
a  digital  system  may  not  be  applicable  in  the  analog  case. 

As  an  alternative  to  fault  simulation,  one  may  adopt  one  of  the  more 
classical  equation  solving  algorithms  for  the  solution  of  1.  Here,  one  first 
measures  m  and  on  the  basis  of  this  measurement,  makes  an  initial  guess  r° 
(usually  taken  to  be  nominal  parameter  vector)  at  the  solution  of  the  equations. 
One  then  evaluates  m°  =  F(r°)  and  compares  it  with  m.  If  m°  =  m,  r°  is  the 
solution  to  the  fault  diagnosis  equation.  If  not,  one  makes  a  new  "educated" 
guess  at  the  solution,  r^ ,  (usually  based  on  the  deviation  between  m  and  m°) 

**Most  industrial  users  of  ATE  obtain  satisfactory  fault  detection  in  digital 
circuits  via  fault  simulation  techniques  but  require  guided  probe  techniques 
in  addition  to  the  fault  dictionary  data  for  fault  diagnosis  (isolation). 
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and  one  repeats  the  process  by  evaluating  =  F(r^)  and  comparing  it  with 
m.  Hopefully,  sequence  of  component  parameter  vectors,  r1 ,  simulated  data 
vectors,  m1  =  F(r1),  is  obtained,  which  "quickly"  converges  to  r  and  m, 
respectively.  Since  the  evaluation  of  F(r1)  is  essentially  equivalent  to  the 
simulation  of  the  system  with  the  faulty  parameter  values,  r1 ,  this  technique 
is  really  another  form  of  fault  simulation.  In  this  case,  however,  one 
simulates  the  system  after  the  data  vector  has  been  measured  and  uses  this 
data  to  make  an  educated  guess  at  a  (hopefully)  small  number  of  parameter 
vectors  at  which  the  system  should  be  simulated.  As  such,  the  approach  has 
been  termed  simulation  after  test  to  distinguish  it  from  the  classical 
approach,  wherein  all  simulation  is  done  before  test. 

At  the  time  of  this  writing,  both  approaches  are  under  study,  neither 
of  which  have  been  shown  to  be  superior.  Fault  "simulation  after  test"  requires 
that  one  include  an  efficient  simulator  in  the  ATE  itself,  which  can  be  used 
for  on-line  computation  of  m1  =  F(r^)  after  the  UUT.  has  been  measured.  On  the 
ohter  hand,  simulation  after  test  eliminates  the  requirement  of  searching  a 
large  fault  dictionary  for  the  (approximate)  data  matches  required  by 
"simulation  before  test".  In  addition,  the  complex  ATPG  requirement  for 
"simulation  before  test"  is  eliminated. 

To  make  "simulation  after  test"  feasible,  however,  an  efficient  equation 
solving  algorithm  is  required  to  obtain  convergence  of  the  r1  sequence  in  a 
reasonable  amount  of  time.  Moreover,  since  "real  world"  failures  in  analog 
circuits  and  systems  often  take  the  form  of  open  and  short  circuited  components 
or  large  parameter  diviations  from  nominal,  the  classical  perturbational 
algorithms  a-la  Newton-Raphson  are  inapplicable.  Fortunately,  in  the  context 
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of  the  fault  diagnosis  problem,  one  can  reasonably  assume  that  relatively 
few  component  parameters  have  failed.  As  such,  even  though  it  is  not  valid 
to  assume  that  r-r°  (the  deviation  of  r  from  nominal)  is  small  in  norm,  it 
is  reasonable  to  assume  that  it  is  small  in  ''rank". 
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Introduction 

After  a  half  century  of  neglect  by  the  electronics  community  the 
past  decade  has  witnessed  an  expanding  effort  in  the  analog  fault  diagnosis 
area.  Indeed,  the  ever  increasing  complexity  of  electronic  circuits  com¬ 
bined  with  the  decreasing  availability  of  trained  maintenance  technicians 
has  pushed  compateji-(Lide.d  testing  (CAT)  to  the  forefront  of  electronics 
research.  Unfortunately,  the  tremendous  strides  which  have  been  made  in 
digital  test  technology  have  not  been  paralleled  by  equal  progress  in  the 
analog  area.  As  such,  even  though  "80%  of  the  boards  are  digital  80%  of 
the  problems  are  analog". 

The  lack  of  progress  in  analog  CAT  vis-a-vis  digital  CAT  may  be 
attributed  to  four  factors: 

i).  the  cost  of  analog  circuit  simulation, 

ii).  the  continuous  nature  of  analog  failure  phenomena, 

iii).  tolerances  on  the  "good"  components  in  analog  circuit, 

iv).  and  the  lack  of  viable  models  for  the  components  in  a  faulty  circuit. 
Moreover,  these  difficulties  have  been  exaggerated  by  the  economics  of  the 
maintenance  environment  which  limits  the  degree  to  which  many  of  the  classical 
tools  of  analog  circuit  design  can  be  used  in  a  CAT  package. 

The  purpose  of  the  present  paper  is  to  describe  a  set  of  criteria  which 
we  believe  a  practical  analog  CAT  algorithm  should  achieve  and  to  indicate 
the  degree  to  which  they  are  met  by  the  various  algorithms  which  have  thus 

*  This  research  supported  in  part  by  the  Joint  Services  Electronic  Program 
at  Texas  Tech  University  under  ONR  Contract  76-C-l 1 36 . 
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far  been  proposed.^  These  criteria  include  computational  requirements, 
numbers  of  test  points  and  test  vectors  employed,  robustness  to  tolerance 
effects,  availability  of  models,  and  the  degree  to  which  the  algorithm  is 
amenable  to  parallel  processing.  Although  many  specific  algorithms  have 
been  proposed  they  may  naturally  be  classified  into  three  categories: 

i).  simulation-before-test, 

ii).  simulation-after-test  with  a  single  test  vector, 

iii).  and  simulation-after-test  with  multiple  test  vectors. 

Each  of  these  three  approaches  to  the  analog  CAT  problem  is  compared  against 
our  criteria,  and,  interestingly,  each  approach  fails  to  meet  at  least  one 
of  the  proposed  criteria. 

Criteria 

A.  Computational  Requirements:  Unlike  a  CAD  algorithm  which  is  used 
only  in  the  initial  design  of  a  circuit  or  system,  a  CAT  algorithm  lives  in 
an  operational  environment  and  thus  must  be  used  repeatedly  each  time  a 
system  fails.  As  such,  a  viable  measure  for  the  computational  cost  of  a  CAT 
algorithm  must  distinguish  between  on-line  computation  which  is  done  in  the 
field  and  must  be  repeated  for  each  unit  under  test  (UUT)  and  off-line  compu¬ 
tation  which  is  independent  of  the  unit  under  test  and  thus  need  only  be  done 
once  at  the  factory  or  a  maintenance  depot.  Indeed,  the  distinction  between 
on-line  and  off-line  computation  is  further  exaggerated  by  high  cost  of 
computing  and  the  dearth  of  trained  personnel  in  a  field  maintenance  environ¬ 
ment  vis-a-vis  that  is  available  at  a  maintenance  depot.  Thus  in  a  CAT  al¬ 
gorithm  a  g Kent  ptiionity  muAt  be.  placed,  on  reducing  the  on-line  computational 
aequOiementA  even  at  the  cost  of  significantly  increasing  the  off-line  com¬ 
putation.  As  such,  an  algorithm  which  is  viable  in  a  design  environment  might 

not  be  acceptable  in  a  maintenance  environment  and  vice-versa.  Indeed,  in  a 
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CAT  algorithm  one  would  be  happy  to  accept  the  cost  of  generating  a  complex 
data  base  in  an  off-line  environment  to  achieve  a  reduction  in  on-line 
computational  requirements. 

B.  Test  Points:  Historically,  analog  circuits  have  been  tested  with 
the  aid  of  a  "bed  of  nails"  tester  which  allows  one  to  make  use  of  test 

data  which  is  not  accessible  via  the  input  and  output  terminals  of  the  circuit 
board.  Unfortunately,  modern  circuit  boards  are  often  multilayered  and/or 
coated,  thereby  limiting  the  applicability  of  the  "bed  of  nails"  concept. 

As  such,  a  modern  CAT  algorithm  must  be  designed  to  work  with  the  test  data 
which  is  available  at  the  externally  accessible  terminals  of  a  printed 
circuit  board.  In  practice,  this  proves  to  be  a  dominating  factor  in  the 
design  of  a  CAT  package,  which  precludes  the  use  of  some  of  the  more 
attractive  algorithms  with  test  point  requirements  which  grow  linearly  with 
circuit  complexity.  In  fact,  circuit  complexity  is  proportional  to  the 
area  of  a  printed  circuit  board  (if  not  a  power  thereof)  while  the  number 
of  accessible  test  points  is  proportional  to  the  edge  length  of  the  board. 

As  such,  in  a  practical  CAT  package  it  is  reasonable  to  require  that  the 
nwmben  of  te*t  point*  gtiow  with  the  *quaAe  Koot  of,  circuit  complexity  (or  less). 

C.  Robustness:  Unlike  a  digital  system  wherein  a  device  is  either  good 
or  bad  in  an  analog  circuit  a  device  is  either  "in-tolerance"  or  "out-of¬ 
tolerance"  and,  as  such,  an  analog  CAT  algorithm  must  be  able  to  cope  with 
the  effects  of  component*  which  one  in-toleAance  but  not  nominal.  Although, 
at  the  time  of  this  writing,  there  is  insufficient  experimental  data  to 
determine  the  import  of  robustness  in  an  analog  CAT  algorithm  it  is,  at 
minimum,  a  factor  of  which  one  must  be  cognizent  and  may,  in  fact,  prove  to 
be  a  dominating  factor  in  the  design  of  a  viable  CAT  package. 
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D.  Models :  Since  most  CAT  algorithms  presuppose  some  form  of  circuit 


simulation  in  their  operation  and  design  of  such  an  algorithm  must  consider 
the  type  and  availability  of  circuit  models  which  are  required  and/or 
available.  In  particular,  does  the  algorithm  use  nominal  circuit  models  or 
6au£te.d  (UacuuX  models?  Indeed,  even  if  nominal  circuit  models  are  used  do 
they  operate  in  their  normal  range?  Finally,  one  must  consider  whether  or 
not  the  algorithm  is  capable  of  dealing  with  "fuzzy"  components  which  do  not 
admit  viable  simulation  models. 

E.  Module  vs.  Parameter  Testing:  Most  analog  fault  diagnosis  algorithms 
can  be  catagorized  as  either  module  oriented  or  parameter  oriented.  In  the 
former  case  the  algorithm  tests  the  input-output  performance  of  the  individual 
modules  or  subsystems  which  make  up  the  UUT  while  in  the  latter  case  the 
algorithm  estimates  a  set  of  parameter  values  which  determine 

the  performance  of  a  given  circuit  component.  Although  one  can  often 
formulate  a  circuit  model  for  a  given  module  thereby  permitting  one  to 
use  a  parameter  oriented  algorithm  to  test  modules,  such  a  process  may 
unnecessarily  complicate  the  test  procedure.  As  such,  a  modulz  ofii&nte.d 
CAT  algofujthn  -ca  pte&Wiid  over  a  parameter  oriented  algorithm  if  it  can  be 
formulated  without  compromising  other  factors. 

F.  In-Si tue  Testing:  Although  secondary  to  the  above  considerations 
the  ideal  CAT  algorithm  should  allow  for  -u t-a-ctue  -tea-tcng.  Since  one  cannot 
control  the  input  signals  applied  to  the  UUT  in-situe  such  an  algorithm 
must  work  with  an  arbitrary  input  signal  rather  than  a  fixed  set  of  test 
vectors. 

G.  Parallel  Processing:  Since  the  CAT  problem  is  inherently  a  large 
scale  systems  problem  it  is  essential  to  exploit  whatever  computational 
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power  is  available  to  reduce  both  on-line  and  off-line  computational 
requirements.  In  particular,  digital  CAT  algorithms  of  ten  use  some  degree 
of  parallel  processing  in  their  implementation.  Given  the  additional 
computational  problems  associated  with  an  analog  CAT  algorithm  the  de.g4.ze 
to  which  an  algorithm  can  be  Implemented  in  parallel  becomeA  a  significant 
factoa  in  determining  its  viability  and  should  therefore  be  included  among 
our  criteria  for  an  analog  CAT  package. 

In  the  above  paragraphs  we  have  described  seven  aspects  of  the  CAT 
problem  which  must  be  considered  in  judging  an  analog  CAT  algorithm. 

Although  we  would  ideally  like  to  formulate  an  algorithm  with  minimal 
computational  requirements  a  moderate  amount  of  off-line  computation  is 
acceptable  since  the  off-line  computation  need  only  be  done  once  and  is 
carried  out  in  a  depot  environment  where  good  computational  facilities  and 
high  level  personnel  are  available.  On  the  other  hand  since  the  on-line 
computation  associated  with  a  CAT  algorithm  is  replicated  for  each  UUT  and 
carried  out  in  a  field  environment  it  must  be  kept  to  a  minimum.  Likewise 
the  test  point  requirements  for  an  analog  CAT  algorithm  must  be  kept  to  a 
minimum.  Although  the  requirement  that  the  number  of  test  points  used  by 
a  CAT  algorithm  grow  with  the  square  root  of  circuit  complexity  is  open  to 
debate  it  is  indicative  of  a  fundamental  limitation  to  the  effect  that  the 
number  of  test  points  should  grow  at  less  than  a  linear  rate  with  circuit 
complexity.  Concerning  the  remaining  criteria  we  want  an  algorithm  that  is 
robust  though  the  significance  of  this  requirement  is  not  fully  understood 
at  this  time.  Similarly,  the  availability  of  circuit  models  to  implement 
an  algorithm  must  be  considered.  Finally,  but  secondary  to  the  above  require¬ 
ments,  it  would  be  desirable  to  have  a  module  oriented  algorithm  which  is 
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amenable  to  in-situe  testing  and  parallel  processing.  These  criteria  are 
sumnarized  in  table  1  along  with  a  set  of  goals  which  one  would  wish  to 
achieve  in  an  "ideal  analog  fault  diagnosis  algorithm". 

CAT  Algorithms 

A.  Simulation-Before-Test:  Although  it  is  essentially  a  brute  force 

search  algorithm  simulation-before-test  is  well  suited  to  the  depot/field 

computational  environment  of  the  CAT  problem  and,  as  such,  it  predominates 

3 

in  most  state-of-the-art  digital  CAT  packages.  On  the  other  hand  its  weak¬ 
nesses  become  more  pronounced  in  the  analog  problem  wherein  it  has  yet  to 
be  successfully  implemented.  Basically,  a  simulation-before-test  algorithm 
is  a  search  algorithm  in  which  one  simulates  the  expected  test  data  which 
would  result  from  various  hypothesized  failures  in  an  off-line  environment. 
Then  when  the  actual  test  data  is  obtained  in  the  field  it  is  compared 
with  the  simulated  results  to  determine  the  failure.  Needless  to  say  the 
technique  requires  immense  amounts  of  off-line  computer  time  to  generate 
the  required  data  base  but  is  extremely  efficient  on-line,  wherein  one  need 
only  compare  the  test  results  with  the  simulated  data  base. 

Unfortunately,  the  cost  of  an  analog  simulation  is  much  greater  than  that 
of  a  digital  simulation.  Moreover,  one  requires  a  much  larger  data  base  in 
the  analog  problem  than  in  the  digital  problem  to  cope  with  the  continuous 
nature  of  the  analog  failure  phenomena  and  the  robustness  problem.  As  such, 
there  is  considerable  doubt  about  the  applicability  of  the  simulation-before¬ 
test  concept  in  an  analog  CAT  package. 

Vis-a-vis  our  criteria  for  analog  fault  diagnosis  simulation-before¬ 
test  requires  extremely  large  amounts  of  off-line  computer  time  but  only  a 
minimum  of  on-line  computer  time.  Additionally,  the  test  point  requirements 
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for  the  algorithm  are  minimal.  On  the  other  hand  the  technique  has  no 
inherent  robustness  and  uses  faulted  simulation  models  for  all  components. 

With  regard  to  the  secondary  factors  the  algorithm  is  module  oriented  and 
amenable  to  parallel  processing  but  not  in-situe  testing.  These  considerations 
are  summarized  in  Table  1. 

B.  Simulation-After-Test  with  a  Single  Test  Vector:  Rather  than  using 
a  search  algorithm  for  fault  diagnosis  one  can  attempt  to  model  the  analog 
fault  diagnosis  problem  as  a  nonlinear  equation  in  which  one  solves  for 
the  internal  variables  or  component  parameters  in  terms  of  the  test  data. 
Although  this  may,  at  first,  seem  to  totally  bypass  the  repetitive  simu¬ 
lation-before-test  algorithm,  a  careful  analysis  will  reveal  that  each  itera¬ 
tion  of  the  required  numerical  equation  solver  amounts  to  a  simulation  of 
the  UUT.  In  this  case,  however,  the  particular  simulations  which  one  carries 
out  are  based  on  known  test  data  rather  than  a-priori  fault  hypotheses.  As 

such,  the  simulations  are  done  on-line  after  the  test  data  has  been  obtained 
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and  the  technique  is  thus  termed  simulation-after-test. 

In  the  case  where  only  a  single  test  vector  is  employed  the  resultant 
fault  diagnosis  equations  are  "almost  linear"  and  may  be  solved  with  the  aid 
of  a  single  (off-line)  sparse  matrix  inversion. The  test  point  require¬ 
ments  for  the  algorithm,  however,  grow  linearly  with  circuit  complexity. 
Interestingly,  this  class  of  algorithms  have  been  discovered  independently 
by  a  number  of  authors  over  the  years ,  most  of  whom  thought  that  they  had 
found  the  "ideal  algorithm"  until  they  fully  appreciated  the  significance 
of  the  test  point  requirement  which  severely  limits  its  applicability. 

From  the  point  of  view  of  our  other  criteria,  however,  the  algorithm  is. 
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indeed,  "ideal".  Off-line  computational  requirements  are  moderate  while 
on-line  computational  requirements  are  minimal.  Moreover,  the  algorithm 
is  inherently  robust  and  requires  no  simulation  models  of  any  kind,  it 
tests  modules,  and  it  is  amenable  to  in-situe  testing.  Finally,  the 
computational  requirements  associated  with  the  algorithm  are  sufficiently 
low  so  as  to  render  the  parallel  processing  question  moot. 

C.  Simulation-After-Test  with  Multiple  Test  Vectors:  One  approach 
to  reducing  the  test  point  requirements  of  the  simulation-after-test 
algorithm  is  to  use  multiple  test  vectors  to  increase  the  number  of 
equations  obtained  from  a  given  set  of  test  points,  thereby  rendering  the 
fault  diagnosis  equation  soluable  with  a  restricted  number  of  test  points. 

The  most  cotimon  form  of  the  multiple  test  vector  algorithm  is  the  multi¬ 
frequency  algorithm  used  in  linear  fault  diagnosis,  though  the  concept 

extends  to  the  nonlinear  case  via  the  use  of  multiple  test  vectors  of  any 

1  2 
type.  ^ 

The  reduced  test  point  requirement  obtained  via  the  use  of  multiple  test 
vectors  is,  however,  achieved  at  the  cost  of  greatly  increasing  the  complex¬ 
ity  of  the  resultant  fault  diagnosis  equations.  Indeed,  the  "almost  linear" 
equations  of  the  single  test  vector  algorithm  are  replaced  by  an  extremely 
complex  set  of  nonlinear  equations  (even  for  linear  systems)  in  the  multiple 
test  vector  algorithm.  Although  these  equations  can  be  made  trackable  in 
the  linear  case  they  appear  to  be  totally  untrackable  in  the  nonlinear  case 
and,  as  such,  most  of  the  advantages  of  the  simulation-after-test  concept 
are  lost  when  multiple  test  vectors  are  employed. 

With  regard  to  our  criteria  the  multiple  test  vector  algorithms  require 
large  amounts  of  on-line  computer  time  though  relatively  little  off-line 
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computer  time  is  required.  In  its  most  obvious  form  the  technique  is 
robust,  though  this  robustness  is  compromised  by  most  of  the  "tricks"  which 
have  been  proposed  to  make  the  multiple  test  vector  fault  diagnosis 
equations  trackable.  Faulted  simulation  models  are  required  and  the  al¬ 
gorithm  is  inherently  parameter  oriented.  Finally,  it  is  not  suited  to 
either  in-situe  testing  or  parallel  implementation. 

Conclusions 

The  above  concepts  are  surmarized  in  Table  1,  wherein  the  various 
criteria,  by  which  an  analog  CAT  algorithm  should  be  measured  are  tabulated, 
the  goals  for  an  ideal  algorithm  are  described,  and  the  degree  to  which  the 
various  algorithms  achieve  these  goals  is  indicated.  From  the  table  it  is 
apparent  that  none  of  the  algorithms  is  fully  acceptable.  Indeed,  even  if 
one  neglects  the  secondary  considerations  regarding  modules  vs.  parameters, 
in-situe  testing,  and  parallel  processing  all  three  approaches  fail  to 
meet  one  or  more  of  the  primary  criteria  (indicated  by  capital  leters  in 
the  table).  As  such,  the  proper  approach  to  the  solution  of  the  analog 
CAT  problem  remains  an  open  question. 
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Table  1:  Performance  of  Analog  Fault  Diagnosis  Algorithms.  Unacceptabl 
performance  factors  are  indicated  by  capital  letters. 
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Analog  Fault  Diagnosis  with  Failure  Bounds* 

C.-c.  Wu,  K.  Nakajima,  C.-L.  Wey,  and  R.  Saeks 
Department  of  Electrical  Engineering 
Texas  Tech  University 
Lubbock,  Texas  79409 


Abstract 

A  simulation-after-test  algorithm  for  the  analog  fault  diagnosis  problem 
is  proposed  in  which  a  bound  on  the  maximum  number  of  simultaneous  failures 
is  used  to  minimize  the  number  of  test  points  required.  The  resultant  al¬ 
gorithm  is  applicable  to  both  linear  and  nonlinear  systems  and  can  be  used 
to  isolate  a  fault  up  to  an  arbitrarily  specified  "replaceable  module". 


*  This  research  supported  in  part  by  the  Joint  Services  Electronic  Program 
of  Texas  Tech  University  under  ONR  Contract  76-C-1136. 
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I.  Introduction 


Conceptually,  analog  fault  diagnosis  algorithms  can  be  subdivided  into 
3 

three  classes;  simulation-before-test,  simulation-after-test  with  a  single 
test  vector,  and  simulation-after-test  with  multiple  test  vectors.  The  former 
is  commonly  employed  in  digital  testing  and  is  characterized  by  minimal  on¬ 
line  computational  requirements.  Unfortunately,  the  high  cost  of  analog 
circuit  simulation  coupled  with  the  large  number  of  potential  fault  modes 
which  must  be  simulated  in  an  analog  circuit  limits  the  applicability  of 
simulation-before-test  algorithms  in  an  analog  test  environment.  As  an 
alternative  to  simulation-before-test,  a  number  of  researchers  have  proposed 
simulation-after-test  algorithms,  in  which  the  internal  system  variables  or 
component  parameters  are  computed  from  the  test  data  via  a  "nonlinear 
equation  solver  -  like"  algorithm.  Indeed,  in  the  case  where  sufficiently 

many  test  points  are  available  only  a  single  test  vector  is  required  and  the 

8  9 

fault  diagnosis  problem  reduces  to  the  solution  of  a  linear  equation.  ’ 

Except  for  the  large  number  of  test  points  required,  this  approach  is  ideally 
suited  to  the  analog  fault  diagnosis  problem  and,  as  such,  a  considerable 
research  effort  has  been  directed  towards  the  problem  of  reducing  its  test 
point  requirements.^  One  such  approach  uses  multiple  test  vectors  to  increase 
the  number  of  equations  obtained  from  a  given  set  of  test  points.  Unfortunate¬ 
ly,  this  is  achieved  at  the  cost  of  greatly  complicating  the  set  of  simultan¬ 
eous  equations  which  must  be  solved  and,  as  such,  the  applicability  of  the 
approach  is  limited. 

The  purpose  of  the  present  paper  is  to  describe  an  alternative  simulation- 
after-test  algorithm  in  which  a  bound  on  the  maximum  number  of  simultaneous 


FILMU) 


t 


35 


FRKHttwo  Pag*  bumc-wot 


failures  is  used  to  reduce  the  test  point  requirements  while  still  retaining 
the  computational  simplicity  inherent  in  a  single  test  vector  algorithm. 

Indeed,  even  though  a  given  circuit  may  contain  several  hundred  components 
it  is  reasonable  to  assume  that  at  most  two  or  three  have  failed  simultaneous¬ 
ly.  As  such,  rather  than  solving  a  set  of  simultaneous  equations  in  n-space 
the  solution  to  our  fault  diagnosis  problem  actually  lies  in  a  two  or  three 
dimensional  submanifold  which  should  yield  a  commensurate  reduction  in  test 
point  requirements.  Unfortunately ,  even  though  we  may  assume  that  at  most 
two  or  three  components  have  failed  we  do  not  know  which  two  or  three,  and  as 
such,  some  type  of  search  is  still  required.  Fortunately,  with  the  aid  of 
an  appropriate  decision  algorithm  the  required  seach  can  be  implemented 
quite  simply. 

Consider  the  circuit  or  system  which  is  illustrated  abstractly  in  figure  1. 


Figure  1.  Test  algorithm  for  abstract  circuit  or  system. 
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Here,  the  individual  circuit  components  or  sub-systems  are  denoted  by 
circles  indexed  from  a  to  n.  These  components  are  subdivided  into  two  groups, 
at  each  step  of  the  test  algorithm,  as  indicated  by  the  dashed  lines  in 
figure  1.  At  each  step  we  assume  that  one  group;  say,  d  through  n; 
is  composed  of  good  components  and  we  use  the  known  characteristics  of 
these  components  together  with  the  test  data  to  determine  whether  or  not  the 
remaining  components;  a,  b  and  c  in  this  case;  are  good.  Of  course,  if 
components  d  through  n  are  actually  good  then  the  resultant  test  results  for 
components  a,  b  and  c  will  be  reliable.  On  the  other  hand,  if  any  one  of  the 
components  d  through  n  is  faulty  the  test  data  on  a,  b  and  c  will  be  unreliable. 
As  such,  we  repeat  the  process  at  the  next  step  of  the  test  algorithm  with  a 
different  subdivision  of  components.  For  instance,  we  may  assume  that  a 
through  d  and  h  through  n  are  good  and  use  their  characteristics  to  test 
components ,  e,  f  and  g.  Finally,  after  a  number  of  such  repetitions  the  test 
results  obtained  at  the  various  steps  are  analyzed  to  determine  the  faulty 
components. 

Of  course,  the  number  of  components  which  may  be  tested  at  any  one  step  is 
dependent  on  the  number  of  test  points  available  while  the  number  of  steps  re¬ 
quired  is  determined  by  the  number  of  components  which  may  be  tested  at  any  one 

step  and  the  bound  on  the  maximum  number  of  simultaneous  failures.  As  such, 
the  procedure  yields  a  natural  set  of  tradeoffs  between  the  numbers  of  test 
points,  simultaneous  failures  and  steps  required  by  the  algorithm.  Indeed, 
since  the  computational  cost  associated  with  each  step  of  the  algorithm  is 

essentially  the  cost  of  a  single  system  simulation  the  latter  parameter  is 

a  natural  measure  of  the  computational  cost. 

In  the  following  section  we  describe  the  simulation  model  used  to  test 
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one  set  of  components  under  the  assumption  that  the  remaining  components 
are  good.  The  model  is  formulated  in  both  the  linear  and  nonlinear  cases 
and  can  be  used  as  readily  to  test  modules  and  subsystems  as  individual 
components.  Moreover,  the  requirement  that  an  appropriate  matrix  be  invertible 
determines  the  maximum  number  of  components  which  can  be  simultaneously  tested 
from  a  given  set  of  test  points  as  well  as  the  allowable  component  sub¬ 
divisions.  In  section  three  two  decision  algorithms  for  analyzing  the  resul¬ 
tant  test  data  are  described.  Indeed  the  required  theory  is  reminiscent, 

though  not  identical  to,  the  t-diagnosibi 1 i ty  theory  developed  for  digital 

4  6 

system  testing  over  the  past  decade.  ’  In  the  context  of  our  application 
we  give  an  exact  decision  algorithm  for  the  case  of  a  single  failure  together 
with  an  analysis  of  the  possible  tradeoffs  between  test  points  and  algorithm 
steps  (read  computer  costs).  Although  an  exact  decision  algorithm  for  the 
multifailure  case  has  yet  to  be  developed  an  heuristic  algorithm  which  is 
applicable  to  both  the  single  and  multifailure  case  is  presented.  The  algorithm, 
which  is  based  on  an  inherently  analog  heuristic,  to  the  effect  that  two 
analog  errors  will  never  cancel,  has  proven  to  be  highly  reliable  while 
simultaneously  reducing  the  number  of  steps  required  from  that  of  the  exact 
(single  failure)  algorithm.  Finally,  section  four  is  devoted  to  a  number  of 
examples.  These  examples  include  linear  circuits  with  12  and  22  components 
which  were  run  on  a  desktop  calculator  and  a  16  bit  mini,  respectively. 

Although  we  have  yet  to  implement  the  algorithm  in  the  nonlinear  case 
the  nonlinear  algorithm  is  identical  to  the  linear  algorithm  except  for  the 
requirement  that  a  nonlinear  simulator  be  used  in  lieu  of  the  linear  simulator. 

As  such,  we  believe  that  the  increase  in  computational  costs  for  the  nonlinear 
case  vis-a-vis  the  linear  case  will  be  in  proportion  th  the  increased  running 
time  for  the  nonlinear  simulator. 
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II .  The  Simulation  Model 

Although  our  test  algorithm  can  be  formulated  in  terms  of  any  of  the 

standard  system  models  for  the  purpose  of  this  exposition  we  will  assume 

2 

a  component  connection  model  for  the  circuit  or  system  under  test.  In  the 
nonlinear  case  the  unit  andeA  teAt  is  represented  by  a  set  of  decoupled  state 
models  characterizing  its  components  and/or  subsystems  together  with  an 
algebraic  connection  equation  as  follows. 


Xj  *  f ( x .  ,a. ) 

1  111  ;  x.(0)  *  0,  1-1,2 . n 


bi  ■  9,(x,.*,) 


(2.1) 


and 


a 


Lnb  +  4iu 


(2.2) 


y  *  l_2-|b  +  L22u  (2.3) 

Here,  a  =  col(a^)  is  the  column  vector  composed  of  the  component  -input  vaAiabZeA, 
b  *  col(b.)  is  the  column  vector  composed  of  component  output  variable*,  u  is 
the  vector  of  exteAnaZ  teAt  input*  applied  to  the  system  and  y  is  the  vector 
of  <su4iem  Adpome*  meaAuned  at  the  vclaLo ua  test  point a.  Although  the  component 
connection  model  is  not  universal  it  is  quite  general  and  subsumes  most  of  the 
classical  topological  connection  models  commonly  used  in  circuit  and  system 
theory.  Moreover,  its  inherently  decoupled  nature  is  ideally  suited  to  the 
test  problem  wherein  we  desire  to  distinguish  between  the  characteristics 
of  the  individual  system  components.  Although  these  components  may  be  taken 
to  be  elementary  RLC  components  and/or  discrete  semiconductor  devices,  in 
practice  the  ''components''  are  taken  to  be  the  "replaceable  modules"  within  the 
circuit  or  system,  under  test;  say,  an  IC  or  a  "throw-away"  circuit  board. 
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At  each  step  of  the  test  algorithm  we  subdivide  the  "components"  into 
two  groups  denoted  by  "1"  and  "2"  with  the  components  in  group  "1"  assumed 
to  be  good  and  used  together  with  the  known  values  of  u  and  y  to  compute  the 
component  input  and  output  variables,  a^  and  b.. ,  for  the  components  in  group 
"2".  Although  computationally  we  prefer  to  work  with  the  decoupled  component 
equations  for  notational  brevity  we  combine  the  equations  for  the  components 
in  each  group  into  a  single  equation 


.fV.aT: 


b1  .gV..1) 


x  (0)  *  0 


(2.4) 


x2  =  f2(x2,a2) 


x*(0)  «  0 


(2.5) 


b2  -  g2(x2,a2) 

Here,  x1,  a1  and  b1  are  the  vectors  of  group  "1"  component  state  variables, 

2  2  2  • 

inputs  and  outputs;  and  similarly  for  x  ,  a  and  b  .  To  retain  notational 
compatibility  with  2.4  and  2.5  we  reorder  and  partition  the  connection 
equations  of  2.2  and  2.3  to  be  conformable  with  2.4  and  2.5  as  follows 


J  a  ,  nhl  .  ,  12.  2  .  .  1  „ 

3  3  L-j  <|  D  ^  L-|  i  D  +  L  -j 


(2.6) 


.2  .  .  21.1  .  ,  22.  2  .  ,  2  „ 
3  -  L^D  +  1  k  +  L^2U 


(2.7) 


y  =  tJ-jb1  +  L^b2  +  L22u 


(2.8) 


Given  equations  2.4  through  2.8  our  goal  is  to  compute  the  group  "2" 

2  2 

component  variables,  a  and  b  ,  given  the  tea*  -input,  u,  the  meoioted 
-'i&jpon.ies ,  y,  and  an  m^umption  to  thz  e^ect  that  thz  gsioup  "I"  ccmponz ti  axe. 
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2  2  -L 

.tot  $awity.  To  this  end  we  assume  that  admits  a  left  inverse,  [L^]  , 

which,  in  turn,  determines  the  allowable  component  subdivisions.  Under 

this  assumption  one  may  then  formulate  a  component  connection  model  for  a 

"paeudo  CsLicllU"  composed  of  the  group  ”1”  components  with  external  input 

vector  up  *  col(u  ,  y)  and  external  output  vector  yp  =  col (a2,  b2)  in  the 

form 


*'  .fV..') 


b1  =  gV,a') 


i  x'(0)  *  0 


(2.9) 


a'  -  K„b’  *  Kj/ 


(2.10) 


■  V'  *  V' 


(2.11) 


Indeed,  some  algebraic  manipulation  of  equations  2.6  through  2.8  together  with 

2  -L 

the  assumption  that  CLgi ^  exists  will  yield 


K11  3 


JW  L.  1 
L11  l11ll21j  l21 


(2.12) 


s12 


i  ^  -  L12a2  '  i12n2  l”'- 

ui2  LnLL2iJ  l22  :  LnLL2iJ 


L21-L22rL2  I'^L 
L11  ‘•11ll21j  l21 


V21 


-0-^ 

LL21J  L21 


(2.13) 


(2.14) 


*22 


I  2  ^22p^2 

l12  l11ll21j  l22 


n2  ]"l*. 
ll21j  l22 


L22[l  ]“*■ 

LnLL2iJ 


ai,rL 


(2.15) 


» 
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Since,  in  our  test  problem  both  u  and  y  are  known,  the  above  equations 

D  2  2 

can  be  solved  via  any  standard  circuit  analysis  code  to  compute  yp  =  a  ,b^). 

Now,  under  our  assumption  that  the  group  "1"  componenta  sre  not  faulty  yp  = 

2  2 

(a  ,b  )  represents  the  inputs  and  outputs  which  actually  appeared  at  the 

terminals  of  the  group  "2"  components  during  the  test.  As  such,  we  may 

determine  which  of  the  group  "2"  components  are  faulty  by  solving  equation 
2 

2.5  with  input  a  and  checking  to  determine  whether  or  not  the  resultant 

2 

output  coincides  with  b  .  Of  course,  since  our  assumption  to  the  effect  that 
the  group  "1"  components  are  not  faulty  may  not  be  valid  -the  neAuLtA  o& 
ttuA  teAt  cuiz  not  x&LLable..  As  such,  we  repeat  the  process  a  number  of 
times  with  different  choices  for  the  subdivision  of  the  components  into 
group  "1"  and  group  "2".  Here,  the  only  constraint  on  the  choice  of  sub- 
divisions  is  the  requirement  that  [L^]  exist  while  the  number  of  com¬ 
binations  employed  is  limited  only  by  the  cost  of  the  required  simulations. 

The  results  of  the  several  steps  in  the  test  algorithm  are  then  analyzed 
via  the  techniques  described  in  the  following  section  to  determine  those 
components  which  are  actually  faulty.  To  this  end  the  results  of  each  step  of 
the  test  algorithms  are  tabulated  as  follows 


Here,  a,b,c,  ...  ,k  denote  the  group  "1  components  for  a  given  step  of  the 
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test  algorithm,  x,y . z  denote  the  corresponding  group  "2"  components 

while  the  binary  annotation  associated  with  the  group  "2"  components  in¬ 
dicates  whether  this  step  of  the  test  algorithm  indicated  that  they  were 
good  (0)  or  bad  (1).  Although  this  tabular  notation  is  somewhat  cumbersome 
we  will  eventually  generate  a  binary  array  indexed  by  the  group  "1"  and 
group  "2"  components  in  the  process  of  our  decision  algorithm  in  which  case 
the  tabular  notation  probes  to  be  convenient. 

For  linear  systems  one  may  formulate  an  identical  algorithm  in  which 
the  component  equations  2.1  are  modeled  in  the  frequency  domain  via 

ai-z1b1  ;  i'1.2,  ---  .n  (2-16) 

where  we  have  suppressed  the  s-variable  for  notational  brevity.  Then  upon 

subdividing  the  components  into  two  groups  characterized  by  the  equalities 
1  a^  1  a^ 

b  =  Z  and  b  =  Z  and  solving  the  resultant  equations  under  the  assumption 
that  O-21 ]~L  exists  one  obtains  an  equation  in  the  form  yp  =  Mup.  Specifically, 

a2  =  M^u  +  (2.17) 


b  =  M„,u  +  M 


21u  M22y 


(2.18) 


where 

u  _  / p|  21  ■  22r,  2  -,-L.  1  -| r i  7 2 f,  11  .  1 2r.  2  -.-L.  1  .  ^  -h2ri_2  i”h  V 

Mn  “  ^hl'hl^hr  L11]{1'Z  !-Lll'L12zL2r  L21;  1  tL21  hrL21J  L22J' 


+  fl2  -l22n2  l”h 

CL12  L11[L21J  L22 


(2.19) 


m  _  /r|  21  . 22  n  2  -t-L.  1  im  7!  ri  11  I  ^2n  2  >”hVh2rL2  l*^l+fL22rL2  1*^ 
m12  (cl1  T ~l1  1  ^21 11  hr  _Z  f-hi-Li2t-L2r  71  1  *-LllLL21  -1  '  ^lrhr 

(2.20) 


M21  =  ‘^L21'i  l’L21{1"zl|-Lll"L121-L21  ^  Ll21}  1  CL21  ~L1 1  (-L21  -1  L22J^"^21]  L 


21  1-nLL-2r  22J '  L  21  22 

(2.21) 
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and 


1P  11  ,  12n  2  -,-L,  1  ,-17lr,  1 2  n  2  ,-L, 


M22  =  -(C41]"‘-L‘in-Z,[.j‘-L^[L‘1]-LL‘1}-|Z'[L^[l21]'L)  -  (2.22) 


Although  these  expressions  appear  to  be  foreboding  they  may  all  be 
computed  with  the  aid  of  only  a  single  matrix  inversion.  Moreover,  since  the 
M.j  are  independent  of  the  test  data  and  computed  in  terms  of  the  nominal 
values  of  the  group  "1"  components  they  may  be  computed  off-line  and  stored  in 
a  data  base  to  be  retreived  at  the  time  a  test  is  conducted.  Furthermore, 
since  only  a  single  test  vector  is  required,  single  frequency  testing  can  be 
employed  in  which  case  the  M. .  need  only  be  computed  at  a  single  frequency. 

'  J 

As  such,  the  only  on-line  computation  required  for  the  fault  diagnosis  of  a 

linear  system  is  the  matrix-vector  multiplication  indicated  by  equations  2.17 

2  2 

and  2.18  together  with  the  computation  of  Z  a  . 

Unlike  the  linear  case,  if  one  is  working  with  a  nonlinear  circuit  or 

2  2 

system,  the  simulations  required  to  compute  a  and  b  require  a-priori  knowledge 
of  y  and  y  and  thus  must  be  carried  out  on-line.  In  practice,  however,  relative¬ 
ly  few  time  steps  are  required  by  these  simulations,  thereby  minimizing  their 
running  time.  Moreover,  all  simulations  are  carried  out  with  nominal  components 
allowing  one  to  use  standard  CAO  circuit  models.  Indeed,  since  the  group  "2" 
component  models  are  only  invoked  at  the  final  step  of  the  analysis  one  can 
avoid  simulating  ''troublesome"  components  by  always  including  them  in  group 
"2"  though  this  usually  means  that  additional  test  points  will  be  required. 

As  such,  one  can  avoid  simulating  "fuzzy”  components  which  do  not  admit  a  viable 
simulation  model  and/or  nonlinear  components.  Indeed,  if  sufficiently  many 
test  points  are  available  to  permit  all  nonlinear  components  to  be  included  in 
group  "2"  a  linear  simulation  model  such  as  that  of  2.17  and  2.18  may  be 


employed  even  for  a  nonlinear  system. 
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III.  Decision  Algorithms 

Since  the  results  of  the  tests  described  in  the  preceding  section  are 
dependent  on  our  assumption  that  the  group  "1"  components  are  not  faulty 
they  are  not  immediately  applicable.  Following  the  philosophy  initiated  by 
Preparata,  Metze,  and  Chein^  in  their  study  of  self  testing  computer  networks, 
however,  if  one  assumes  a  bound  on  the  maximum  number  of  faulty  components 
it  is  possible  to  determine  the  actual  fault(s)  from  an  analysis  of  the  test 
results  obtained  at  the  various  steps  in  the  algorithm.  To  this  end  we  will 
give  a  complete  analysis  of  the  theory  required  to  locate  a  single  fault 
together  with  an  heuristic  which  is  applicable  to  the  multiple  fault  case. 

Let  us  assume  that  at  most  one  circuit  component  is  faulty  and  that  the 
test  results  obtained  from  a  given  step  of  the  algorithm  indicate  that  all 
group  "2"  components  are  good  as  indicated  in  the  following  table. 


In  this  case  we  claim  that  the  group  "2"  components  are,  in  fact,  good.  Indeed, 
if  a  group  two  component  were  actually  faulty  then  our  test  results  are  in¬ 
correct,  which  could  only  happen  if  one  of  the  group  "1"  components  was 
faulty.  As  such,  the  system  would  have  two  faulty  components  contradicting 
our  assumption  to  the  effect  that  at  most  one  component  is  faulty. 

Now,  consider  the  case  where  the  results  from  a  given  step  of  the  test 
algorithm  indicate  that  exactly  one  group  "2"  component  is  faulty;  say,  x 
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In  this  case  the  same  argument  we  used  above  will  guarantee  that  the  components 
which  test  good;  say,  y  through  z;  are,  in  fact,  good.  On  the  other  hand  we 
have  no  information  about  x.  It  may  be  faulty  or,  alternatively,  the  test 
result  may  be  due  to  a  faulty  group  "1"  component. 

Finally,  consider  the  case  where  two  or  more  group  "2"  components  test 
bad  in  a  given  step  as  indicated  in  the  following  table. 


Since,  under  our  assumption  of  a  single  failure,  it  is  impossible  for  two  or 
more  group  "2"  components  to  be  faulty,  this  test  result  implies  that  at  least 
one  of  the  group  "1"  components  is  bad.  On  the  other  hand  since  we  have  assumed 
that  there  is  at  most  one  faulty  component  the  faulty  group  "1”  component  is 
the  only  faulty  component  and,  as  such,  the  group  "2"  components  are  all  good. 

Consistent  with  the  above,  at  each  step  of  the  test  algorithm,  either  all 
or  all  but  one  of  the  group  "1"  components  are  found  to  be  good.  As  such,  if  we 
choose  our  subdivisions  so  that  the  components  which  are  founa  to  be  good  at 
one  step  of  the  algorithm  are  included  in  group  "1"  in  all  succeeding  steps 
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we  will  eventually  arrive  at  a  group  "l",  all  of  whose  components  are  known 
to  be  good.  As  such,  the  test  results  obtained  at  that  step  will  be  reliable, 
thereby  allowing  us  to  accurately  determine  the  faulty  components  in  group  '‘2". 
Although  the  number  of  components  in  group  “1"  and  group  "2"  may  vary  from  step 
to  step  (especially  if  we  work  with  multivariate  components)  if  we  assume  that 
group  "1"  contains  n-m  components  n-m  components  and  group  "2"  contains  "m" 
components  at  each  step  of  the  algorithm  then  the  process  will  terminate  in 
approximately  n/m  steps.  Since  the  computational  cost  of  the  algorithm  is 
proportional  to  the  number  of  steps  (essentially  the  cost  of  one  simulation 
per  step)  while  m  is  determined  by  the  number  of  allowable  test  points  the 
ratio  n/m  represents  a  natural  measure  of  the  possible  tradeoffs  between  test 
point  and  computer  requirements  when  employing  the  algorithm  in  a  single  fault 
made. 

Unlike  the  single  fault  case,  at  the  time  of  this  writing,  we  do  not  yet 
have  an  exact  decision  algorithm  for  the  multiple  fault  case.  Following  Liu, 
however,  theproblem  can  be  greatly  simplified  if  one  adopts  an  "analog 
heuristic"  to  the  effect  that  two  independent  analog  failures  will  never 
cancel.5  Needless  to  say,  this  is  an  inherently  analog  heuristic  since  two 
binary  failures  have  a  fifty-fifty  chance  of  canceling  one  another.  In  the 
analog  case,  however,  two  independent  failures  are  highly  unlikely  to  cancel 
one  and  another  (as  long  as  one  works  with  reasonably  small  tolerances). 

Recall  from  our  discussion  of  the  single  fault  case  that  whenever  a  test 
result  indicates  that  a  component  is  good  then  it  is,  in  fact,  good.  Although 
this  is  not  rigorously  true  in  the  multiple  failure  case  it  is  true  under  the 


47 


assumption  of  our  heuristic.  For  instance,  consider  the  test  results  indicated 
in  the  following  table  in  which  x  is  found  to  be  good. 


H 


i 


i 

i 


:  t 


Now,  if  x  is  actually  faulty  there  must  be  a  faulty  group  "1"  component  whose 
effect  is  to  cancel  the  error  in  x  as  observed  during  this  step  of  the  test 
algorithm.  This  is,  however,  forbidden  by  our  heuristic  and,  as  such,  we 
conclude  that  x  is  actually  good. 

Interestingly,  our  heuristic  can  be  carried  a  step  further  than  indicated 
above  since,  under  our  heuristic,  a  bad  group  "1"  component  would  normally 
yield  erroneous  test  results.  An  exception  would,  however,  occur  if  some  of 
the  group  "1"  components  are  totally  decoupled  from  some  of  the  group  "2" 
components.  As  such,  if  prior  to  our  test  we  generate  a  coupling  table  (by 
simulation  or  a  sensitivity  analysis)  which  indicates  whether  or  not  a  faulty 
group  "1"  component  will  effect  the  test  results  on  a  group  "2"  component,  our 
heuristic  may  be  used  to  verify  that  certain  group  "1"  components  are  good 
whenever  a  good  group  "2"  component  is  located.  Consider  for  example  the 
following  table 
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in  which  a  "1"  in  the  i-j  position  indicates  that  the  test  results  for  component 
i  are  affected  by  component  j  while  a  "0"  in  the  i-j  position  indicates  that 
component  j  does  not  affect  the  test  results  for  component  i.  Now,  since 
component  x  has  been  found  to  be  good  in  this  test  our  heuristic  implies  that 
those  group  "1"  components  which  are  coupled  to  x  in  this  test  shows  that  z  is 
good  the  heuristic  implies  that  b  and  c  are  also  good.  This,  with  a  single 
test  we  have  verified  that  x,z,a,b,c,  and  k  are  all  good. 

Since  in  any  practical  circuit  the  coupling  table  is  composed  mostly  of 
l's  it  has  been  our  experience  that  relatively  few  steps  of  the  algorithm  will 
yield  a  complete  diagnosis.  To  implement  the  heuristic,  however,  one  must 
assume  that  the  maximum  number  of  faulty  components  is  strictly  less  than 
the  number  of  group  two  components.  If  not,  the  test  results  at  each  step 
may  show  that  all  group  "2"  components  are  faulty,  in  which  case  no  reliable 
test  information  is  obtained.  Moreover,  the  degree  to  which  the  number  of 
group  "2"  components  exceeds  the  maximum  number  of  faulty  components  determines 
the  number  of  algorithm  steps  which  will  be  required  to  fully  diagnose  a 
circuit. 

Although  no  exact  decision  algorithm  for  the  multiple  failure  case 
presently  exists  it  is  noteworthy  that  the  underlying  combinatorial  decision 
problem  is  quite  similar  to  the  t-diagnosibility  problem  usually  associated  with 
self  testing  computer  networks,  wherein  the  multiple  fault  problem  has  been 
resolved. In  that  problem,  however,  one  computer  tests  another  with  the 
rest  of  the  network  being  decoupled,  whereas  in  our  problem  a  subset  of 
components  test  all  remaining  components  since  there  exists  no  practical 
mechanism  for  decoupling  components  in  an  analog  circuit  or  system.  In  any 
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event  the  problem  ere  similar  and.  In  fact,  W.'  baa  already  Wated 
a  generalization  of  the  t-diagnosibility  problem  in  which  one  aubaet  of 
computers  on  a  network  tests  another.  Us  such,  we  believe  that  an  exact 
decision  algorithm  for  our  problem  can  be  formulated  in  the  multiple  failure 

case. 


i 


IV.  Examples 


To  illustrate  the  exact  decision  algorithm  for  the  single  fault  case 
consider  a  system  composed  of  eight  components;  a,b,  ...  ,h;  in  which  any 
five  may  test  the  remaining  three.  Initially,  we  let  a  through  e  represent 
the  group  "1"  components  and  f,g  and  h  represent  the  group  "2"  components 
and  assume  that  the  test  results  for  this  first  step  are  as  indicated  in 
the  following  table. 


Employing  our  exact  algorithm  for  the  single  fault  case  the  above  table  in¬ 
dicates  that  components  f  and  g  are  good  and,  as  such,  we  move  them  into 
group  "1"  for  the  second  step  of  the  algorithm  obtaining 


Since  this  test  indicates  that  two  group  "2"  components  are  bad  which  contradicts 
our  single  fault  assumption  the  faulty  component  must  be  in  group  "1"  implying 
that  d,e,  and  h  are  all  good.  We  therefore  move  these  components  into  group  "1" 
and  implement  the  final  step  of  our  algorithm  in  the  form 


\ 

ii  2** 
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c 
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Since  all  group  "1"  components  are  known  to  be  good  this  final  test  is 

reliable  and  indicative  of  the  fact  that  b  is  the  faulty  component. 

2 

Note,  that  the  requirement  that  be  left  invertible  may  make  it  im¬ 
possible  to  use  some  component  subdivisions  in  which  case  an  alternative 
sequence  of  steps  may  be  required  in  the  above  process.  For  instance,  if 
h,e,d,f,  and  g  is  not  an  allowable  subdivision  the  last  step  in  the  above 
process  might  be  replaced  by 


indicating  that  c  and  h  are  good.  Now,  a  final  step  in  which  c,e,d,f,  and  g 
make  up  group  "1"  will  be  reliable  as  indicated  below. 


Now,  consider  the  same  single  fault  example  in  which  our  heuristic  algorithm 
is  applied  using  the  coupling  table  indicated  below 
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According  to  our  heuristic  f  and  g  are  good  and,  moreover,  everything  in  group 
"1"  which  is  coupled  to  either  f  or  g  is  good.  As  such,  we  conclude  from  this 
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first  step  that  f,g,a,c,d,  and  e  are  all  good.  Thus,  taking  group  "1"  to  be 
e,d,f,g  and  a  in  the  next  step  will  yield  a  reliable  test  for  b,d  and  h  as 
above. 

Finally,  consider  the  case  where  at  most  two  failures  are  assumed  with 
the  first  step  in  our  test  algorithm  yielding: 
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Consistent  with  cur  heuristic  f,a,d,  and  e  are  found  to  be  good  in  this  step. 
Incorporating  these  components  into  group  "1"  for  the  following  step  we  obtain: 
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which  gives  us  no  information  in  the  multiple  failure  case*.  As  such,  we  try 
another  allowable  combination  obtaining  the  following  table 
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indicating  that  h.f.g.a,  and  e  are  good.  Coupled  with  our  previous  knowledge 


Actually  one  can  deduce  that  at  least  one  of  the  group  "1"  components  is 
bad  since  all  three  group  "2"  components  cannot  be  bad  via  our  two  fault 
assumption.  This,  in  turn,  implies  that  at  most  one  of  the  group  "2" 
components  is  bad. 
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hat  d  is  good  this  implies  that  all  group  "l"  components  are  good  and  hence 
this  last  step  in  our  algorithm  reliably  indicates  that  b  and  c  are  the 
faulty  components. 

To  obtain  more  realistic  examples  the  above  techniques  are  applied  to 
the  12  and  22  component  linear  amplifier  circuits  shown  in  figure  2  using 
simulated  test  data  for  various  numbers  of  simultaneous  failures,  choices 
of  test  point  locations,  and  both  decision  algorithms.  All  analysis  for  the 
12  component  circuit  was  done  on  an  HP  9825  desktop  calculator  while  the  22 
component  examples  were  run  on  a  TI  990/20  minicomputer.  The  results  of 
some  150  simulations  of  the  algorithm  are  tabulated  in  table  1.  where  the 
number  of  test  points,  simultaneous  faults,  and  the  decision  algorithm  employed 
are  indicated.  The  results  of  the  various  simulations  are  indicated  by  the 
ambiguity  of  the  resultant  diagnosis.  For  instance,  in  our  simulation  of  the 
12  component  circuit  with  3  test  points,  one  failure  and  the  exact  algorithm 


Figure  2.  a)  12  component  amplifier  and  b)  22  component  amplifier. 

All  stages  of  the  amplifier  circuits  have  nominal  op-amp 
gains  of  1.6,  nominal  resistances  of  10K  ohms,  and 
nominal  capacitance  values  of  .001 yf  while  the  feedback 
capacitors  have  nominal  values  of  lOOpf. 


12  runs  were  made  (one  with  each  component  faulty).  On  10  occasions  the 
fault  was  located  exactly  while  the  fault  was  located  exactly  while  the  fault 
was  located  up  to  an  ambiguity  set  composed  of  two  components  on  2  occasions. 
Finally,  we  note  that  the  5th  run  of  the  12  component  circuit  indicated  by  an 
asterisk  in  the  table  represents  a  simulation  in  which  the  good  components 
were  set  at  +/- 2%  off  of  nominal  to  test  the  robustness  of  the  algorithm. 
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V.  Conclusions 


Although  the  proposed  algorithm  is  still  new  and  we  are  just  beginning  to 
investigate  its  performance  in  the  nonlinear  case  the  algorithm  ptomises  to 
meet  most  of  the  criteria  set  in  reference  7.  Although  the  on-tint  computa¬ 
tional  nzquOi&mzntA  for  the  algorithm  do  not  compare  with  a  simulati on-before¬ 
test  algorithm  they  can  be  kept  within  reasonable  bounds.  Indeed,  unlike 
most  simulation-after-test  algorithms  no  iterative  on-line  computation  is 
required.  Moreover,  one  can  limit  the  on-line  computation  by  restricting  the 
number  of  algorithm  steps  (at  the  price  of  increasing  the  ambiguity  in  the 
resultant  diagnosis).  Furthermore,  in  the  linear  case  and/or  in  the  case  where 
there  are  sufficiently  many  test  points  available  to  permit  all  nonlinear  and 
fuzzy  components  to  be  included  in  group  "2"  the  major  part  of  the  computation 
required  by  the  algorithm  can  be  done  off-line. 

In  general  the  proposed  algorithm  permits  one  to  tradeoff  between  on-line 
computational  requirements  and  tzit  point*.  Indeed,  as  indicated  in  table  1., 
one  can  reduce  the  tzii  point  azquOizmcnti  to  quite  reasonable  levels  though 
this  is  usually  achieved  at  the  cost  of  increasing  the  number  of  steps  in  the 
algorithm  (and  hence  its  on-line  computational  requirements).  In  particular, 
our  simulations  indicate  that  the  algorithm  comes  close  to  achieving  the  /n 
test  point  goal  set  in  reference  7. 

With  respect  to  the  remainder  of  the  criteria  specified  in  reference  7 
the  algorithm  "looks  good".  In  particular,  all  simulations  are  carried  out 
using  nominal  component  model*,  it  can  test  linza/i  and  nonlinzan  module*  of 
arbitrary  size,  and  is  amenable  to  in-iituz  testing  and  panallzl  pnocz*J>ing 
techniques  (since  several  steps  of  the  algorithm  can  be  carried  out  simul¬ 
taneously)  . 


At  the  present  time  the  major  open  question  with  respect  to  the  perform¬ 
ance  of  the  algorithm  is  its  Tcfaoit:>tes4 .  Indeed,  there  is  nothin  in  the 
algorithm  to  make  it  inherently  robust  though  our  initial  test  for  robustness 
indicated  by  the  asterisk  in  table  1.  proved  to  be  favorable. 
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Analog  Fault  Simulation 

C.-c.  Wu**,  A.  Sangiovani-Vincentelli*.  and  R.  Saeks** 

I .  Introduction 

After  a  half  centry  of  neglect  by  the  circuits  and  systems  conmunity  the 
past  decade  has  witnessed  the  emergence  of  a  research  effort  in  the  analog 
circuit  maintenance  area.  The  various  algorithms  which  have  been  thus  far 
proposed  for  the  analog  fault  diagnosis  problem  may  naturally  be  subdivided 
into  two  classes  termed  "simulation-before-test"  and  "simulation-after-test". 
The  former  are  commonly  used  in  digital  system  test  algorithms  and  require 
an  automatic  test  program  generator  (ATPG)  which  simulates  the  responses  of 
"all  possible"  failures.  This  is  typically  done  at  a  maintenance  depot  with 
the  simulated  responses  being  recorded  and  shipped  to  the  field  where  the 
response  of  the  unit  under  test  (UUT)  is  compared  with  the  simulated  responses 
to  determine  the  failure.  The  major  advantage  of  simulation-before-test  is 
that  it  is  ideally  matched  to  the  depot/field  maintenance  environment  with  the 
largest  part  of  the  computation  done  only  once.  As  such,  the  technique  is 
ideally  suited  for  digital  testing  where  the  binary  nature  of  the  problem 
keeps  the  number  of  failures  to  be  simulated  within  bounds  and  eliminates 
tolerence  problems.  Unfortunately,  in  the  analog  problem  we  must  cope  with 
a  continuum  of  possible  failures  and  simultaneously  deal  with  good  components 
which  are  in  tolerance  but  not  nominal.  As  such,  a  tremendous  number  of  simu¬ 
lations  are  required  by  a  simulation-before-test  algorithm,  while  some  type 
of  decision  algorithm  is  required  to  cope  with  the  tolerance  effects. 

*  Dept,  of  Elec.  Engrg.  and  Comp.  Science,  Univ.  of  California  at  Berkeley, 
Berkeley,  CA  90024. 

**  Dept,  of  Elec.  Engrg.,  Texas  Tech  Univ.,  Lubbock,  TX  79409.  This  research 
supported  in  part  by  the  Joint  Services  Electronics  Program  at  Texas  Tech 
University  under  ONR  Contract  76-C-1136. 
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Unlike  simulation-before-test,  simulation-after-test  uses  an  "equation 
solver-like"  algorithm  to  compute  the  parameters  of  the  UUT  components  in  the 
field.  Since  most  such  algorithms  require  iterative  evaluation  of  the  equation 
co  be  solved,  the  UUT  is  effectively  simulated  at  each  iteration,  though  the 
simulation  is  based  on  actual  test  data  rather  than  hypothesized  failure  data. 
The  simulation  process  is,  thus,  carried  out  after  testing  the  UUT  and  hence 
the  choice  of  terminology.  The  advantage  to  such  an  approach  is  that  the 
faulty  component  parameters  are  computed  explicitly,  thereby,  eliminating  the 
ambiguity  caused  by  the  use  of  discrete  simulation-before-test  data  and 
tolerance  effects.  Although  relatively  few  simulations  are  required  for  each 
UUT,  they  must  be  carried  out  in  the  field  rather  than  the  depot  and  they  must 
be  repeated  for  each  UUT. 

The  purpose  of  the  present  paper  is  to  describe  a  research  effort  directed 
at  alleviating  some  of  the  difficulties  in  developing  a  simulation-before-test 
algorithm  for  analog  fault  diagnosis.  The  underlying  philosophy  and  motivation 
for  our  formulation  is  discussed  in  section  2,  along  with  a  derivation  of  the 
required  differential-interpolative  fault  diagnosis  formula.  Finally,  section  3 
is  devoted  to  a  number  of  illustrative  examples  of  the  approach.  These  include 
both  linear  and  nonlinear  examples  formulated  in  the  frequency  and  time  domains, 
respectively. 
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1 1 .  A  Differential-Interpolati ve  Algorithm 

Although  any  practical  fault  diagnosis  algorithm  must  be  able  to  handle 
systems  with  a  hundred  or  more  components,  from  an  intuitive  point  of  view 
our  algorithm  is  best  illustrated  in  the  two  component  cases  where  the 
parameter  space  can  be  displayed  graphically.  Say,  we  are  dealing  with  an 
RC  circuit  for  which  the  parameter  space  is  illustrated  in  figure  1. 


Figure  1:  Parameter  space  for  RC  circuit. 


Here,  R  and  (1  represent  normalized  parameter  values,  wherein,  the  nominal 
parameter  values  are  transformed  to  the  origin.  In  the  most  general  simulation- 
before-test  algorithms  one  assumes  that  the  faulty  parameter  values  may  lie 
anywhere  in  the  R-C  plane  and  therefore  carries  out  simulations  along  a  two 
dimensional  discrete  array  spread  over  the  entire  plane. 

Fortunately,  in  a  "real  world"  testing  environment  one  can  assume  that 
only  a  "limited  number  of  components"  fail  simultaneously.  In  our  two  component 
example  we  may  therefore  assume  that  either  R  or  C  has  failed  with  the  other 
remaining  nominal  in  which  case  the  circuit  need  only  be  simulated  at  a  discrete 
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set  of  points  along  the  coordinate  axes  in  the  R-C_  plane  denoted  by  x's  in 
figure  1.  As  such,  the  number  of  simulations  required  is  significantly  de¬ 
creased.  Indeed,  this  is  one  of  the  major  advantages  of  the  simulation-before¬ 
test  concept  as  compared  to  simulation-after-test  algorithms  which  typically 
fail  to  exploit  a  "limited  number  of  failures"  assumption. 

While  the  above  described  approach  has  been  used  with  considerable  success 
in  digital  system  testing,  wherein,  the  axes  are  binary  and  no  tolerance  problems 
are  encountered,  it  is  not  well  suited  to  the  analog  test  problem.  First,  an 
analog  failure  may  occur  anywhere  along  the  axis  and  hence  some  type  of  approx¬ 
imation  scheme  is  required  to  interpolate  between  the  discrete  simulations. 
Secondly,  a  "good"  component  is  assumed  to  be  in-tolerance  though  it  may  not  be 
nominal.  As  such,  in  an  analog  environment  the  "limited  number  of  failures" 
assumption  implies  that  the  solution  to  our  fauK  diagnosis  problem  lies  near, 
but  not  necessarily  on,  the  coordinate  axes  as  indicated  by  the  shaded  regions 


Figure  2:  a)  Solution  space  under  a  single  failure  assumption, 
b)  Illustration  of  the  differential-interpolative 
diagnosis  algorithm. 

While  we  might  choose  to  simply  fill  the  shaded  region  with  additional 
simulations,  the  cost  of  such  a  process  may  prove  to  be  excessive.  Rather,  we 
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exploit  the  fact  that  the  deviations  of  good  component  parameters  from 
nominal  are  small  and  use  a  1st  order  Taylor  series  approximation  to  approximate 
the  deviation.  We  note  that  such  an  approach  cannot  be  used  to  locate  the  faulty 
parameter  values  which  may  be  far  from  nominal;  indeed,  it  is  often  infinite 
or  zero;  though  it  can  be  used  to  cope  with  the  tolerance  effects. 

Our  differential-interpolative  approach  thus  uses  a  classical  minimum  distance 
algorithim  to  locate  the  general  region  of  the  faulty  parameter  values  indicated  by 
the  circle  in  figure  2a  (which  is  magnified  in  figure  2b).  Now,  it  is  assumed  that 
the  simulated  values  of  the  system  responses;  f -j ,  f,,,  and  f3;  corresponding  to 
the  points;  c-j ,  c ,,,  and  c3;  are  available  along  with  the  associated  inverse 
sensitivity  matrices;  J^1,  j^1 ,  and  J'1 .  We  then  interpolate  these  data  points 
to  approximate  the  system  responses  and  the  associated  inverse  sensitivity 
matrice  along  the  axis  by  functions  f(c)  and  J(c)"1.  Although  any  interpolation 
can  be  employed  we  have  had  our  best  results  using  a  bilinear  interpolation 
for  f  (which  gives  exact  results  in  the  linear  case) 'and  a  second  order  polynomic 
interpolation  of  J”  .  Now,  if  x  denotes  the  faulty  parameter  vector  and  m  de¬ 
notes  the  measured  system  responses  then  a  1st  order  Taylor  series  approximation 
combined  with  our  interpolation  will  yield  the  (approximate)  equality 

m  =  f(c)  +  J(c)[x  -  c]  1 . 

for  those  values  of  c  near  x.  Equivalently, 

[x  -  c]  =  J(c)-1[m  -  f ( c ) ]  2. 

Interestingly,  by  invoking  the  Projection  theorem  one  can  reduce  the  above 
vector  equation  to  a  scalar  equation  and  simultaneously  eliminate  the  requirement 
for  storing  the  inverse  sensitivity  matrices.  Indeed,  the  vector  [x-c]  will  be 
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perpendicular  to  the  axis  at  the  point  c  which  makes  the  closest  approach  to 
the  fault.  As  such,  if  ec  denotes  the  unit  vector  in  the  direction  of  the 
axis  then 

0  =  e*[x  -  c]  =  e*J(c)"1[m  -  f(c)]  3. 

which  can  be  solved  for  the  faulty  parameter  value,  c.  Note,  our  goal  is  to 
solve  for  c,  not  x,  since  we  are  interested  in  locating  the  faulty  parameter 
value  in  the  presence  of  the  tolerance  problem,  but  we  really  do  not  care  to 
compute  the  deviations  from  nominal  in  the  good  parameters. 

To  summarize,  if  rather  than  simply  storing  the  simulated  circuit  responses, 
f . ,  we  also  store  the  vectors  e£jT^  then  the  tolerance  effects  associated  with 
the  good  components  can  be  completely  removed  from  our  fault  diagnosis 
algorithm  -  at  least  up  to  the  approximation  error  induced  by  the  interpolation 
process  and  Taylor  series  expansion.  Since  most  good  circuit  simulation  codes 
include  a  package  for  generating  sensitivity  matrices  at  little  additional 
cost  over  and  above  that  involved  in  simply  simulating  the  circuit  responses 
the  approach  can  be  implemented  with  only  a  minimal  increase  in  simulation 
costs.  As  such,  the  major  expense  associated  with  the  approach  lies  with  the 
storage  requirements  (for  the  f .  and  e^JT^  vectors)  which  are  approximately 
double  that  of  a  classical  fault  simulation  algorithm. 

Although  the  above  derivation  has  been  illustrated  in  the  two  dimensional 
case  with  a  single  faulty  parameter  it  can  be  readily  extended  to  a  general 
setting,  say  with  several  hundred  components  and  three  or  four  simultaneous  faults. 
If  one  assumes  p  simultaneous  faults  then  p  inner  products  are  required  to  apply 
the  Projection  theorem  yielding  p  equations  and  p  unknowns  to  be  solved  for 
the  faulty  parameter  values.  Otherwise  the  formulation  for  the  general  case  is 
identical  to  the  single  fault  case  described  above. 
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III.  Examples 


In  this  section,  three  examples  are  given,  two  of  them  for  linear  systems 
and  one  for  the  nonlinear  case.  All  of  these  examples  were  simulated  on 
an  HP9825A  programmable  calculator,  and  yielded  fairly  good  results. 

Our  first  example  is  a  second  order  low  pass  filter.  The  filter  contains 
five  components,  K,  ,  and  Cj,  while,  the  circuit  diagram  is  shown  in 

figure  3-1 
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(3-4) 


0\s)  = 


•k[S4c1c2r1  +  SC  +  sc2j 


—  (r  s)  = 


-K[S2C1R2R1  +  SR^l-K)] 


(3-5) 


6f  _  -KCS  C, Rt Ro  +  SR,  +  SR,] 

^  (r,5)  -  - M2 - 2 - J —  (3-6) 

where  D  =  S2C1C2R1R2  +  S[R2C2  +  R]C2  +  R^O-K)]  +  1 

Since  we  have  five  parameters  in  the  transfer  function,  five  distinct  test 
frequencies  are  required  to  provide  sufficient  information  for  diagnosis.  The 
f(r)  vector  is  then  formed  by  the  transfer  functions  corresponding  to  each 
frequency. 


f(r)  = 


f(r,s1) 

f(r,S2) 

f(r,S3) 

f(r,S4) 

f(r,S5) 


(3-7) 


Similarly,  the  sensitivity  matrix  is  given  by 


6k  ^R^’V 


'  6C^r’V 


J(r)  = 


6K 


(3-8) 


•  &2^V 
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Next,  we  generate  a  data  base  in  which  f(r)  and  erJ(r)  are  evaluated  at 
selected  points  on  each  axis.  Since  the  bilinear  interpolation  gives  exact 
results  in  the  linear  case,  we  need  not  choose  too  many  points  on  each  axis, 
in  fact,  for  this  simulation  example  15  points  were  selected  along  each  axis. 

Given  the  measured  data  for  a  faulty  circuit,  the  processor  searches  the 
data  base,  to  locate  the  point  a2  which  minimizes  ||m-f(a)  j] .  We  then  choose 
the  two  points  a-j ,  on  both  sides  of  a2  along  the  same  axis,  and  use  these 
three  points  to  evaluate  the  coefficients  of  the  bilinear  polynomial 
in  the  interval  (c^,  a^)  also,  we  use  e^J'^a^),  e^J'^a^)  and  e^J~^(a3)  to 
compute  the  coefficients  of  the  second  order  polynomial  approximate  to  e£<J’^(a) 
Finally,  a  golden  section  linear  search  is  used  to  solve 

e£rT(a)(f(a)-m)  =  0  (3-9) 

The  faulty  diagnosis  results  are  listed  in  table  3-1.  Here,  the  nominal 
values  of  K,  ,  R2>  Cj,  C2  are  1.6,  IKn,  IKfi,  0.16yF  and  0.16yF  respectively  and 
the  faulty  parameter  is  underlined  in  the  table 


Table  3-1 


1 

2 

3 

5 

6 

K 

0.6 

1.62 

1.62 

1.58 

1.62 

4_ 

R1 

1090 

2500 

1090 

1070 

1050 

1050 

R, 

930 

1040 

10 

930 

930 

1045 

C 

0. 163y 

0 . 1 61  yi 

0.157y 

0.25y 

0.1 57y 

0 . 162y 

0. 162y 

0. 162y 

0. 162y 

0.157y 

0.25y 

0 . 162y 

Result 

K 

R1 

R2 

C1 

C2 

K 

0.591 

2492 

19.4 

0.239y 

4.01 
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In  the  first  simulation,  K  is  the  faulty  component  with  a  value  of  0.6, 
the  other  four  components  are  5%  or  so  off  their  nominal  values,  the  simulation 
result  shows  that  K  failed,  and  locates  it  at  K=0 . 591 .  The  same  remarks  apply 


to  the  other  five  simulations. 

Although  the  technique  generally  yields  satisfactory  results  occasional 
errors  occur  when  the  good  components  are  too  far  out  of  tolerance.  For  instance, 
the  following  parameter  values  K- T . 62 »  R^  =  1070,  R2=910,  C-|=0.5w,  and  C2=0.172y 
led  to  an  erroneous  result.  The  simulation  shows  that  C2  is  failed  with  the 
value  of  0.179p.  However,  the  faulty  component,  in  this  simulation,  is 
actually  C-j.  If  we  sketch  a  two  dimensional  representation  of  the  C-| ,  C2  plane 
the  difficulty  becomes  clear.  Figure  3-2  shows  that  C2  is  too  far  away  from 
its  own  nominal  value,  and  thus  instead  of  locating  the  error  at  a  we  expect 
the  simulation  result  locates  the  failure  at  B,  with  the  differential  term  still 
pointing  toward  the  actual  failure  denoted  by  x. 


Figure  3-2 

Our  second  example  is  a  fourth  order  low  pass  filter,  comprised  of  two 
cascade  second  order  low  pass  filters,  the  circuit  diagram  is  shown  in  Figure  3-3. 
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Figure  3-3 


The  procedure  for  generating  the  data  base  is  similar  to  the  one  described 
above,  except  here,  ten  different  frequencies  are  required.  The  fault  diagnosis 
results  are  described  in  table  3-2.  Here,  the  nominal  parameter  values  are 
Kj  =  1-6,  R]  =  R2  ■  IK,  C1  =  C2  =  0.16u,  K2  =  1.2,  R3  =  R4  =  1  - 5K,  C3  = 

-  0.2u,  while  the  faulty  parameter  is  again  underlined  in  the  table. 

Our  final  nonlinear  example  is  composed  of  a  diode  loaded  by  a  shunt  RC 
circuit  as  illustrated  in  Figure  3-4.  The  diode  is  modeled  by  the  character¬ 


istic  function 


(3-10) 


Figure  3-4 

Now,  instead  of  working  with  frequency  domain  transfer  functions,  we  work 
in  the  time  domain.  A  state  equation  for  this  circuit  is  given  by 


V 


(V.-VJ/V 


t  n  c 

-  '  /  -  UF" 


(3-11) 
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Table  3-2 


1 

2 

3 

4 

5 

6 

K1 

1.62 

1.62 

1.57 

1.57 

0.36 

1.61 

R1 

1070 

1070 

1055 

1050 

1070 

1043 

r2 

920 

920 

1030 

300 

930 

1010 

C1 

0.162m 

0.157m 

0.157m 

0.162m 

0.161m 

0.11m 

C2 

0.161m 

0.162m 

0.162m 

0.157m 

0.158m 

0.157m 

k2 

1.17 

3,45 

1.22 

1.22 

1.21 

1.18 

R3 

600 

1550 

1560 

1430 

1430 

1535 

R4 

1560 

1570 

1450 

1570 

1590 

1478 

C3 

0.21m 

0.21u 

0.195m 

0.205m 

0.195m 

0.202m 

C4 

0.205m 

0.195m 

0.12u 

0.203m 

0.203m 

0.198m 

Result 

R3 

K2 

C4 

R2 

K, 

i 

C1 

3.546 

0.125m 

339.7 

0.357 

0.118 
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The  goa1  is  to  integrate  this  differential  equation  so  as  to  build  a 
V  (r)  vector  and  f(r)  vectors  as  in  the  previous  examples. 

Numerical  techniques  can  be  used  to  compute  V  ( t )  at  any  instant  t. 

In  this  example,  the  V  (r)  vector  was  elevated  by  applying  the  fourth  order 
Runge-Kutta  method.  Note  that  since  there  are  four  independent  parameters, 

R,  C,  IQ,  and  V-p  in  equation  (3-11)  V  (r,  t)  should  be  evaluated  at  four 
different  time  instants  to  build  a  V  (r)  vector. 

The  sensitivity  matrix  is  generated  via  the  linear  differential  equations 


*Vc 

II 

fD 

< 

1 

»,)/vt  ev 

0 

0 

6Vt 

I_ 

=  (V 

VjC  c 

%  (Vi-vc)/vT 

-  Vi)e  1  c  1  +  E 

6V 

V 

6V 

c 

w 

-  -j-+  E 

R  u 

c 

5R 

(3-12) 

(3-13) 

(3-14) 


(3-15) 


Where  E  = 


V/VT 


J_ 

RC 


6V.  . 

Again,  a  fourth  order  Runge-Kutta  method  is  applied  to  obtain  -^(r,  t), 
«VC  .  6Vc  .  6Vc  .  6I0 

and  5C~^r’  at  the  f°ur  specified  instants. 

The  simulation  results  are  suirmerized  in  table  3-3.  The  nominal  values  of 


IQ»  Vj,  R,  C  are  0.2,  0.1,  IK,  and  0.25  respectively. 
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Table  3-3 


1 

2 

3 

4 

5 

!o 

0.035 

0.21 

0.2 

0.22 

0.002 

VT 

0.11 

5.3 

0.11 

0.11 

0.09 

R 

920 

930 

6250 

1070 

1090 

C 

0.23 

0.27 

0.255 

2.4 

0.23 

I. 

K 

R 

C 

I 

Result 

0 

0 

0.0312 

5.449 

6001 

2.678 

0.003 
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IV.  Conclusions 

For  the  simulation-before- test  approach  to  fault  diagnosis,  we  gain 
from  the  fact  that  most  computation  can  be  done  by  off-line  computation,  thus 
greatly  reducing  the  repetitive  on-line  computation  associated  with  many 
fault  diagnosis  algorithms.  From  a  practical  point  of  view,  the  economics  of 
such  an  approach  are  extremely  attractive.  Unfortunately,  the  simulation-before¬ 
test  approach  is  subject  to  a  certain  degree  of  ambiguity  introduced  by  good 
components  which  are  in-tolerance  but  not  nominal. 

In  this  paper,  we  have  proposed  a  simulation-before-test  algorithm  for 
analog  fault  diagnosis,  in  which  a  differential-interpolative  technique  is 
used  to  eliminate  the  ambiguity  caused  by  tolerance  effects.  Our  approach 
has  been  tested  with  satisfactory  results  in  both  the  linear  and  nonlinear 
cases.  In  fact,  for  the  linear  case,  the  approach  provides  an  exact  inter¬ 
polation  for  f(c)  on  the  axes,  and  thus  reduces  the  amount  of  simulation-before¬ 
test  data  required  on  each  axis.  Although  this  is  not  true  for  nonlinear  case, 
the  diagnosis  results  are  still  very  attractive.  Of  course,  occasional  errors 
may  occur  when  the  good  components  are  too  far  out  of  tolerance.  This  phenomena 
is,  however,  expected  and  well  understood.  Indeed,  the  difficulty  occurs  only 
when  the  1st  order  Taylor  series  approximation  is  too  good.  Because  this 
phenomena  will  rarely  occur  in  the  real  world,  we  believe  that  it  may  be  neg¬ 
lected  in  a  practical  algorithm. 
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Multitest  Diagnosibil i ty  of  Nonlinear  Circuits  and  Systems 
A.  Sangiovanni-Vincentell i*  and  R.  Saeks** 

I.  Introduction 

During  the  past  decade  a  considerable  research  effort  has  been  devoted 
to  the  analog  fault  diagnosis  problem  wherein  one  desires  to  locate  faulty 
circuit  components  given  the  overall  circuit  response  to  one  or  more  test 
vectors.  Conceptually  the  process  may  be  described  by  a  nonlinear 
equation 

y  -  f (a ,u )  1. 

where  y  represents  the  measured  response  to  the  test  vector  u  given  the 
faulty  parameter  vector,  a.  Since  u  is  known  and  y  is  a  measureable 
quantity  the  fault  diagnosis  problem  may  be  resolved  by  simply 
solving  1.  for  a  given  u  and  y.  Unfortunately,  in  practice,  the  dimension 
of  y  is  limited  by  the  number  of  accessible  test  points  in  the  circuit 
and  is  typically  smaller  than  the  dimension  of  the  parameter  vector, 
thereby  precluding  the  direct  solution  of  1 .  To  alleviate  this  difficulty 
a  set  of  test  vectors;  {u^u^,  ...  ,unl;  is  employed  yielding  the  set  of 
simultaneous  equations 

yi  =  f(a,u . )  ;  i  =  l  ,2,  ...  ,m  2. 

Since  the  parameter  vector,  a,  is  independent  of  the  choice  of  test  vector 
this  process  effectively  increases  the  number  of  available  equations  without 

*Dept.  of  Elec.  Engrg .  and  Comp.  Sci.,  Univ.  of  California  at  Berkeley, 
Berkeley,  Ca.  90024. 

**Dept.  of  Elec.  Engrg.,  Texas  Tech  Univ.,  Lubbock,  TX.  79409.  This  research 
supported  in  part  by  the  Joint  Services  Electronics  Program  at  Texas  Tech 
Univ.  under  ONR  contract  76-C-l 1 36 . 
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increasing  the  number  of  unknowns.  More  Concisely,  if  we  let  =  col(y^) 
and  f  (a)  =  col (f (a, u^))  the  "multitest"  fault  diagnosis  problem 
reduces  to  the  solution  of 

1  -  F  (a)  3. 

Needless  to  say  once  equation  3.  has  been  formulated  its  solution  is 
amenable  to  standard  algorithms.  The  problem,  however,  is  to  determine 
whether  or  not  there  exists  a  set  of  test  vectors  (u-i^,  ...  .u^}  such  that 
3.  is  solvable  in  an  appropriate  sense.  To  this  end  we  will  formulate  a 
diagnosibi 1 ity  criterion  directly  in  terms  of  the  function  f  which  determines 
the  degree  to  which  the  equation  y.  =  F  (a)  will  be  solvable,  given  an 
"optimal"  choice  of  test  vectors.  Since  this  criterion  is  a  property  of 
the  circuit  rather  than  the  test  algorithm  it  can  therefore  be  used  as  a 
design  aid  with  which  to  choose  test  points  and/or  to  aid  in  designing 
"testable  circuits". 

,  3 

The  authors  have  previously  formulated  such  criteria  for  linear  circuits 
and  for  memoryless  nonlinear  circuits,^  respectively.  The  present  formulation 
for  dynamical  nonlinear  circuits  is  obtained  by  simply  reformulating  the 
memoryless  nonlinear  theory  in  an  infinite  dimensional  setting.  Surprisingly, 
however,  since  the  parameter  vector,  r,  remains  finite  dimensional  the 
infinite  theory  yields  a  cUm&nt>-LonaZ  £&i>t  mcutALx.. 

Let  U  denote  a  (topological)  space  of  admissible  test  inputs  and  let  Y 
denote  a  Hilbert  space  of  test  outputs.  We  then  model  a  general  nonlinear 
dynamical  system  as  a  function  f:UxA  -»  Y  where  y  =  f(u,a)  denotes  the  response 
of  the  system  to  an  input  u  given  the  parameter  vector  a  e  A.  Here,  A  is  an 

t/ 

open  subset  of  R  which  defines  our  parameter  space  and  f  is  assumed  to  be 


SO 


t 


n 

continuous  in  u  and  continuously  (Frechet)  differentiable  in  a .  The  problem 

at  hand  is  to  solve  for  a  given  a  set  of  input/output  pairs  (u-,y^), 

i=l,2 . m.  We  say  that  the  system  is  locally  dLa.gnoia.ble.  at  a  point  oQ  c  A 

if  there  exists  an  open  neighborhood,  V,  of  aQ  such  that  for  every  a  e  V 

f(u,a)  f  f (u ,aQ)  for  some  u  e  U.  Similarly,  we  say  that  the  system  is 

locally  dlagno&able.  if  it  is  locally  diagnosable  for  almost  all  aQ  e  A  (in 

the  sense  of  Lebesque  measure).  Finally,  if  M(a)  is  a  matrix  valued  function 

of  a  e  A  we  say  that  a  point  aQ  is  a  acgulaa  point  if  there  exists  an  open 

neighborhood,  V,  of  aQ  such  that  M(a)  has  constant  rank  in  V. 

af 

To  define  our  test  matrix  we  let  J(u,o)  =  (u,a)  denote  the  Frechet 

oCl 

2 

derivative  of  f  with  respect  to  a  evaluated  at  a  fixed  u  and  o.  With  u  e  U 


fixed,  f  maps  A  to  Y  and  hence  J(u,a)  is  a  linear  transformation  mapping  from 
k  2 

R  to  Y  .  Moreover,  for  k  =  1  such  a  mapping  may  be  represented  by  an  element 
a  ^ 

of  Y  hence  --  (u,a)  may  be  identified  with  an  element  of  Y  and 
3ai 


J(u,a)  =  row  (u ,a) ]  e  Yk 


Now  let  w  denote  a  positive  measure  defined  on  the  Borel  sets  of  U  such  that 
w(v)  >  0  for  every  open  set  V  contained  in  U.  Finally,  we  introduce  the  test 
matrix 


Rw(a)  =  /  J*(u,a)J(u,a)dw(u)  5. 

L 

Recall  that  J(u,a)  is  a  mapping  for  R  to  Y  hence  its  adjoint  is  a  mapping 

L 

from  Y  to  R  and  thus  Rw(a)  is  a  matrix  valued  function  even  though  U  and  Y 

2 

may  be  infinite  dimensional.  Also  note  that  the  self  duality  of  Y  is  necessary 
for  Rw(a)  to  be  well  defined. 


i 
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Theorem  1 :  Let  w  be  any  admissible  measure  for  which  Rw(a)  exists 
for  all  a  c  A.  Then  the  system  is  locally  diagnosable  at  a  regular 
point  aQ  e  A  if  and  only  if  R^a)  is  nonsingular. 

Proof:  By  invoking  the  integral  form  of  the  mean  value  theorem  we  have 

1 

f(u,a)  -  f(u,a  )  =  /  J (u ,ta  +  (l-t)a)dt  [a  -  an]  6. 

0  o  0  0 

for  any  a  in  a  neighborhood  of  aQ  and  u  e  U.  Now,  assume  that  the  system  is 
not  locally  diagnosable  at  aQ  in  which  case  there  exists  a  sequence  {a^ }  e  A 
approaching  ag  for  which 

ffu.aj)  =  f(u,aQ)  7. 

Letting  a  =  in  6.  yields 

1 

0  =  /  0(u,tao  +  (l-t)a.j  )dt  la.  -  a  J  8. 

1 

=  /  0(u,ta  +  (l-t)ajdt  [a.]  *  9. 

0  u  ii 

where  a.  =  [a.  -  a„]/||[a.  -  a  ]||.  Since  a.  is  normalized  to  lie  in  the  unit 
1  1  0  11  1  0  11  1 
U 

sphere  of  R  which  is  a  compact  set  a.  admits  a  convergent  subsequence  a. 

1  \ 

2 

whose  limit,  a,  also  has  unit  norm.  Using  the  convergent  subsequence  in  9. 
we  have 

0  =  U(u,ta  +  (l-t)a.j  )dt  [a.  ]  =  /  J(u,B.)dt  [a]  =  0(u,on)a  10. 

0  0  'k  1k  0  0 

Since  a  has  unit  norm  it  is  non-zero  while 

a*R(«Ja  =  /  a*J*(u,a_)0(u,o)a  dw(u)  =  0  11. 

W  O  jj  o 
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implying  that  R  (o  )  is  singular  since  it  is  a  positive  semidefinite  matrix 

W  0 

2 

which  is  not  positive  definite. 

Conversely,  if  Rw(a0 )  is  singular  it  follows  from  the  assumption  that 
aQ  is  a  regular  point  and  lemma  1  of  reference  1  that  there  exists  an  open 
neighborhood  V  of  aQ  and  a  continuous  R  valued  function  c(a)  /  0  defined  on 
V  such  that 

0  =  C*(a)R  (a)c(a)  =  /  C*(a)J*(u,a)J(u,a)c(a)  dw(u)  12. 

W  U 

Since  w(v)  >  0  and  J(u,a)  is  continuous  in  a  this  implies  that 

0(u,o)c(a)  =  0  a  e  V  13. 

Finally,  we  define  a  curve  a(s)  c  V  by  the  differential  equation 

=  c{a)  ;  a(0)  =  aQ  14. 

Substituting  a(s)  into  f(u,a)  and  computing  its  derivative  with  respect  to 

2 

s  via  the  chain  rule  we  obtain 

~{u,a(s) )  =  J(u,a(s))|f  =  J(u,a(s))c(a)  =  0  15. 

showing  that  f(u,a(s))  is  constant  along  a  curve  emanating  from  oQ.  Since 
a(s)  is  independent  of  u  this  implies  that  the  system  is  not  locally 
diagnosable  at  aQ  thereby  completing  the  proof.  □ 

Note,  the  proof  uses  the  finite  dimensionality  of  the  parameter  space  but 
does  not  require  that  U  and  Y  be  finite  dimensional.  As  such,  the  theory 
is  equally  valid  for  memoryless  and  dynamical  systems.  Furthermore,  the 
result  is  independent  of  the  choice  of  the  measure  w. 
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To  use  the  theorem  as  a  test  for  local  diagnosability  we  may  invoke 
Theorem  2  of  reference  1.  Alternatively,  one  may  make  the  mild  assumption 
that  f (u,a)  is  (real)  analytic  in  a  in  which  case  the  hypotheses  of  that 
theorem  are  automatically  satisfied  and  the  following  result  is  obtained. 

Theorem  2.  Assume  that  f(u,a)  is  analytic  in  a  and  let  w  be  any  ad- 

missable  measure  for  which  Rw(a)  exists  for  all  a  e  A.  Then  the  system 

is  locally  diagnosable  if  and  only  if  R  (a  )  is  nonsingular  for  some 

w  o 

°o  E  A* 

Proof:  Let 


p  =  max  (rank[R  (a)]}  16. 

aeA  W 

which  must  be  achieved,  say  at  a,  since  rank[Rw(a)]  takes  on  only  finitely 

many  values.  Now,  let  M(a)  be  a  p  by  p  submatrix  of  Rw(a)  which  has  rank 

p  at  a  and  consider  det[M(a)].  Since  rank[M(a)]  =  p  det[M(a)]  f  0  verifying 

that  det[M(a) J  is  not  identically  zero.  This  function  is,  however,  analytic 

in  a  (since  f(u,a)  is  analytic  in  a)  and  since  it  is  not  identically  zero, 

2 

its  zero  set  is  nowhere  dense  in  A.  Letting 

BC  =  {a  e  A,  det[M(a)]  =  0}  17. 

denote  the  closure  of  the  zero  set  and  B  denote  the  complement  of  Bc  the 

hypotheses  of  theorem  2  of  reference  1  are  satisfied.  Indeed,  since  Bc  is 

2 

the  closure  of  a  nowhere  dense  set  it  has  Lebesque  measure  zero  while  the 
fact  that  det  M[(a)]  f  0  for  a  e  B  implies  that 


p  >_  rank  [Rw(a)3  >_  rank[M(a)]  =  p  ;  a  e  B 


18. 


hence 


i 


« 


rank  [Rw(a)3  =  p  ;  a  e  B 


19. 


Since  the  hypotheses  of  Theorem  2  of  reference  1  are  satisfied  in  our  case 
the  conclusion  of  the  theorem  to  the  effect  that  the  system  is  locally 
diagnosable  if  and  only  if  Rw(aQ)  is  nonsingular  for  some  ao  follows. 

Consistant  with  the  proof  of  Theorem  2  rank[Rw(a)]  is  a  generic  property 
(i.e.,  constant  almost  everywhere)  and  hence  we  may  refer  to  p  as  the  generic, 
tank  of  Rw(«)  when  f(u,a)  is  analytic  in  a.  As  such,  one  may  check  the  non¬ 
singularity  of  Rw(“)  with  a  randomly  chosen  aQ  and  the  result  of  the  theorem 
may  be  restated  via: 

Corollary  1:  Assume  that  f(u,a)  is  analytic  in  a  and  let  w  be  any 

admissible  measure  for  which  R  (a)  exists  for  all  a  e  A.  Then  the 

W 

system  is  locally  diagnosable  if  and  only  if  the  generic  rank  of 

Rw(a)  is  k. 

Once  it  has  been  verified  that  a  system  is  locally  diagnosable  it  still 

remains  to  pick  a  set  of  test  signals  {u.|,u2 . uml  and  to  solve  the 

resultant  set  of  equations 


JL  = 


y^  f  (u-j  ,a) 
y2  f(u2,a) 


F(a) 


f<V> 


20. 


for  a  e  A.  For  this  purpose  we  require  that 
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-  f£<O0> 


J(u1  ,aQ) 

¥,<V“o> 

-  £k(V“o> 

J(u2,aQ) 

. 

ll  <u2*“o> 

*  *  * 

Sk(“n,-o> 

21 


have  column  rank  k  in  which  case  a  Newton-Rapson  type  of  algorithm  will  be 
assured  of  converging  to  aQ  from  a  sufficiently  good  initial  guess.  More 
generally,  if  this  Jacobian  matrix  has  a  generic  column  rank  of  k  such  an 


algorithm  will  converge  to  almost  any  solution  of  17.  given  a  sufficiently 

3f 

good  initial  guess.  Recall  from  equation  1.  that  each  ^  (u^.c^)  is  an 
element  of  Y  and  hence  we  are  dealing  with  an  m  by  k  matrix  whose  entries 
are  elements  of  Y. 

Theorem  3:  Let  w  be  an  admissible  measure  for  which  R  ,(a)  exists  for 

'  W 

all  a  e  A  and  assume  that  Rw(aQ)  is  nonsingular  at  a  regular  point 

aQ  e  A.  Then  there  exists  a  sequence  u.  e  1);  i=l  ,2,  ...  ,  nw  k; 

for  which  (a  )  has  column  rank  k. 

9a  0 

Proof:  If  |—  (u,a  )  =  0  for  all  u  then  the  first  column  J(u,a  1  is  zero  for 
9a-|  0  0 

all  u  implying  that  the  first  row  and  the  first  column  of  R  (a  )  is  zero. 

W  U 

This,  however,  contradicts  the  assumption  that  R  (an)  is  nonsingular  and  we 

w  u 

3  F 

may  therefore  assume  that  there  exists  a  u,  e  U  for  which  (u.,a  )  f  0.  As 

I  C7  CX 1  1  O 

such,  there  exists  u-j  e  U  such  that  has  a  column  rank  of  at  least  1.  Using 
this  fact  as  the  starting  point  for  an  inductive  hypothesis  we  assume  that 
there  exists  u.;  i  =  1,2,  ...  n;  such  that  the  matrix  Hn  has  column  rank 
p  <  k  where  n  p.  We  now  desire  to  verify  the  existence  of  a  vector  un+-j  e  U 
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for  which  the  corresponding  matrix  H  .j  has  column  rank  greater  than  or 
equal  to  p  +  1 .  To  this  end  we  let  T  be  a  nonsingular  matrix  of  scalars 

which  operate  on  the  columns  of  Hn  in  such  a  way  that  the  (p  +  l)st  column 

through  the  kth  column  of  HnT  is  zero.  Since  T  is  nonsingular  will 
have  column  rank  greater  than,  or  equal  to,  p  +  1  for  some  un+-j  if  and  only 
if  Hn+.jT  has  column  rank  greater  than,  or  equal  to,  p  +  1.  Because  of  the 
special  form  of  HnT,  however,  this  will  be  the  case  if,  and  only  if,  the 
bottom  row  of  Hn+1T  given  by 

fa(Vr°o)T  *  J(“n+l’“o)T  22' 

is  non-zero  in  columns  p  +  1  through  k  for  some  u^ .  If  this  is  not  the 

case  we  may  let  t  denote  the  (p+l)st  column  of  T  which  is  non-zero  since  T 

is  nonsingular,  in  which  case  we  have 

J(u,aQ)t  =  0  23. 

for  all  u  e  U.  This,  however,  implies  that 

t*Rj«Jt  =  /  t*J*(ua  )J(ua  )t  dw(u)  =  0  24. 

w  o  u  u  u 

which  contradicts  the  assumption  that  ^(“q)  is  nonsingular.  As  such,  there 

must  exist  a  un+-|  for  which  Hn+1  has  column  rank  greater  than,  or  equal  to 

3F 

p  +  1.  Repeating  the  argument  inductively  until  an  Hm  =  |^(a0)  with  column 
rank  k  is  obtained  now  completes  the  proof  of  the  theorem.  Note  that  since 
n  <  p  at  each  step.m  <_  k. 

Our  purpose  in  the  above  has  been  to  indicate  a  mechanism  by  which  the 
existing  diagnosibility  theory  for  linear  systems^  and  memoryless  nonlinear 
systems  can  be  extended  to  the  general  nonlinear  dynamical  case.  Indeed,  the 
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I 


same  tools  can  be  used  to  extend  the  theory  in  any  number  of  directions.  For 

instance,  one  might  use  the  rank  of  as  a  measure  of  the  degree  to  which 

a  circuit  fails  to  be  diagnosable.  Indeed,  the  resultant  mzaMine.  oh 

tzAtabiJLLty  is  a  natural  generalization  of  the  frequency  domain  measure  of 

testability  introduced  in  reference  4.  Moreover,  the  frequency  domain 

5 

criterion  for  choosing  test  signals  introduced  by  Priester  and  Clary  may 
also  be  extended  to  the  case  of  general  nonlinear  systems  via  the  above 
described  formulation.  Finally,  a  vector  space  criterion  for  diagnosability 
similar  to  that  of  reference  3.  may  be  formulated. 
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I.  Introduction 

Historically,  symbolic  network  analysis  has  been  motivated  by  the  prob¬ 
lems  of  circuit  design  and,  as  such,  the  emphasis  has  been  placed  on 

quickly  and  efficiently  obtaining  a  symbolic  transfer  function  from  a 

2  3 

given  set  of  circuit  specifications.  ’  In  an  operational  or  maintenance 

environment,  however,  one  is  typically  given  a  prescribed  nominal  circuit 

and  desires  determine  the  effect  of  various  (possibly  large)  perturbations 

thereon.  This  is  the  case  in  a  power  system  where  one  is  given  a  fixed 

network  and  desires  to  determine  the  effect  of  proposed  modifications  thereto. 

Alternatively,  in  the  problem  of  analog  circuit  fault  diagnosis  one  desires 

to  simulate  the  effect  of  a  number  of  alternative  failures  to  compare  the 

a 

simulated  data  with  the  observed  failure  data. 

In  such  an  operational  or  maintenance  environment  numerous  perturbations 
of  tne  nominal  circuit  are  studied  and,  as  such,  significant  computational 
efficiencies  can  be  obtained  if  one  first  generates  a  data  base  in  temps 
of  the  nomina1  circuit  parameters  and  then  extracts  the  appropriate  symbo’ic 
transfer  function  from  tne  data  base  each  time  a  different  symbolic  transfer 
is  reauired.  Of  course  the  benefit  to  be  achieved  via  such  an  approach  is 
dependent  on  the  size  of  the  data  base  and  the  ease  with  which  a  symbolic 
transfer  function  may  be  retrieved  therefrom. 

♦This  research  supported  in  part  by  the  Joint  Services  Electronics  Program 
at  Texas  Tech  University  under  QNR  Contract  76-C-l 136. 
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The  obvious  manner  in  which  to  generate  such  a  data  base  is  to  simply 
pre-compute  the  coefficients  of  all  required  symbolic  transfer  functions 
and  store  them  in  the  data  base.  Retreival  from  such  a  data  base  is,  of 
course,  immediate  but  the  data  base  may  become  overly  large.  Indeed,  the 
number  of  transfer  functions  which  must  be  stored  is  0(k^-)  where  k  is  the 
total  number  of  potentially  variable  circuit  parameters  and  £  is  the  maximum 
number  of  circuit  parameters  which  may  vary  simultaneously.  An  alternative 
approach  is  to  store  the  nominal  transfer  function  information  and  then  use 
Householder's  formula^  to  compute  the  required  symbolic  transfer  functions. 
In  such  a  data  base  we  need  only  store  Q(n  )  transfer  functions  where  £  is 
the  total  number  of  component  output  terminals  but  retreival  requires 
0(n3+p3)  multiplications  where  p  is  the  actual  number  of  circuit  parameters 
wnich  vary  simultaneously.  Since,  in  practice,  n  >>  p  the  retreival  process 
reauires  approximately  O(n^)  multiplications  and  is  dominated  by  the  large 
dimensional  matrix  mul tipi ication  reouired  by  Householder's  formula  rather 
than  the  lew  dimensional  inverse. 

In  the  present  paper  we  will  formulate  an  alternative  data  base  for 

2 

the  symbolic  transfer  functions  which  also  requires  0(n  )  entries,  Dut  for 
which  retreival  requires  only  O(p^)  multiplications.  Since  o  is  typically 
small  this  is  tantamount  to  immediate  retreival. 

In  the  remainder  of  this  introduction  we  will  review  the  properties  of 
the  component  connection  model  for  a  large  scale  circuit  or  system^  which 
serves  as  the  starting  point  for  our  theory.  The  data  base  and  retreival 
formulae  for  the  case  where  p  <_  2  are  formulated  in  section  2.  while  the 
general  retreival  formula  is  derived  in  section  3.  Section  4.  is  devoted 


94 


to  the  problem  of  retreiving  sensitivity  formulae  from  the  data  base  while 
section  5.  deals  with  the  problem  of  updating  the  data  base  when  the  nominal 
circuit  parameters  are  changed.  Finally,  section  6.  is  devoted  to  an  example 
illustrating  the  theory. 

The  component  connection  model  is  an  algebraic  model  for  an  inter¬ 
connected  dynamical  system  which  subsumes  the  classical  topological  models 
but  is  more  readily  manipulated  both  analytically  and  computationally.  The 
motivation  and  justification  of  the  model  are  discussed  in  detail  in 
reference  1  and  will  not  be  repeated  here.  The  component  connection  model 
takes  the  form  of  the  set  of  simultaneous  equations 

b  *  Z(jja  1 .1 


and 


=  Lnb  ♦  L 1 2U 


l2lb  ♦  L22u 


1.2 

1.3 


Here,  Z  (=Z(jj))  is  a  frequency  dependent  matrix  characterizing  the  decoupled 
system  components  with  composite  component  input  and  output  vectors  a  and  b, 
respecti vely .  On  the  other  hand  the  L . i ;  i , j  *  1,2;  matrices  are  frequency 
indeDendent  connection  matrices  characterizing  the  coupling  between  the 
comoos ite  component  vectors,  a  and  b,  and  the  comDosite  system  input  and 
output  vectors,  u  and  y,  respectively . 

A  little  algebra  with  the  component  connection  equations  will  readily 
reveal  that 


S  =  L 


22 


L21 ( 1  ‘ 


ZLll)' 


'ZL 


12 


1.4 
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.  1 


where  S  ( =  S(joa))  is  the  composite  system  transfer  function  matrix 
characterizing  the  external  behavior  of  the  system  via 

y  =  S(jui)u 


1.5 


Often,  rather  than  working  with  the  entire  S  matrix  we  find  it  convenient 
to  work  with  its  individual  entries;  sqv,  q  *  1,2,  ...  ,£  and  v  =  1,2,  . . . , v : 
which  are  related  to  the  component  connection  model  via 


,qv  ,  ,qv  ,q  _  ,,  x-lz,v 
s  i22  +  L21U  )  4l12 


Here  l^2  is  the  q-v  entry  in  1_22;  q  =  1,2,  ...  ,£  and  v  =  1,2,  . . . ,  v ;  is 

the  qth  row  of  ;  q  =  1,2,  ...  ,£;  and  L^2  is  the  vth  column  of  L12; 
v  =  1,2,  ...  ,  v. 

Finally,  since  we  are  interested  in  analyzing  the  effects  of  perturbing 
one  or  more  components  from  their  nominal  values,  we  decompose  Z  into 
nominal  and  perturbation  terns  in  the  form 


1.6 


q 


where 


Z  ’  Zo  +  Z1 


7  _  r  k.k  k 

Z,  *  /.  c  Sr 

1  k*l 


1 .7 


1.3 


Here,  c  (=  c  (j->))  is  a  column  vector,  r  (*  r  (jj)  is  a  row  vector,  and 


.  K 


is  the  scalar  perturbation  for  the  kth  potentially  variable  component 


U  |/  L 

oarameter.  In  a  typical  application  one  is  given  c  ,  r  ,  and  5  ; 

k  =  1,2,  ...  ,]<;  characterizing  k_  potentially  variable  component  parameters 
though  at  most  £  such  parameters  vary  in  any  given  analysis;  P  i  £  1  k. 

Indeed,  p  <<  k_  in  most  applications.  Finally,  we  note  that  Z-j  can  be  expressed 
more  concisely  in  the  form 
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where 


Z  =  CaR 


1  '  2' 
C  =  [c1  ,  cc. 


. :  Cp  ] 


The  above  described  notation  formulated  for  the  component  connection  model 


is  summarized  in  table  1. 
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latri  x 

Type 

Dimension 

Index 

a 

composite  component  input 
vector 

m  x  1 

- 

b 

composite  component  output 
vector 

£  x  1 

- 

u 

composite  system  input  vector 

£  x  1 

- 

y 

composite  system  output  vector 

£  X  1 

- 

Ln 

connection  matrix 

m  x  n 

- 

l21 

connection  matrix 

£  X  £ 

- 

L21 

qth  row  of  l-2i 

1  x  £  q  = 

1,2,  ...  ,£ 

l12 

connection  matrix 

m  x  v_ 

- 

l12 

vth  column  of 

m  x  1  v  = 

1,2,  ...  ,v 

L22 

connection  matrix 

1  *  1 

- 

q-v  entry  in 

lxl  q  =1 ,2, 

...  ,£i  v=l ,2,  . . 

s 

composite  system  transfer 
function  matrix 

a  x  i 

• 

sqv 

q-v  entry  in  S 

1  x  1  q=l,2. 

. . .  ,£;  v  =  l ,2, 

Z 

composite  component  transfer 

£  x  m 

- 

Zo 

nominal  composite  component 

n  x  m 

- 

Z1 

composite  component  transfer 
function  perturbation  matrix 

n  x  m 

k 

c 

column  vector  characterizing 
perturbation  of  kth  parameter 

n  x  1  k  = 

1,2,  ...  ,k. 

C 

1/ 

array  of  the  c  vectors  for  the 
parameters  chich  actual  vary 
(  row[ckl) 

£  X  P 

k 

r 

row  vector  characterizing 

1  x  n  k  = 

1,2,  ...  ,k 

R 

array  of  rk  vectors  for  the 
parameters  which  actually 
vary  (col[rk]) 

p  x  m 

6k 

kth  variable  parameter 

lxl  k  * 

1,2,  ...  ,k 

a 

array  of  5k,s  for  parameters  P  x  p 

which  actually  vary  (diag[6Kl) 

Table  1.  Summary  of  Component  Connection  Model 
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II.  The  Data  Base 


Our  data  base  is  composed  of  the  following  family  of  (frequency 
dependent)  scalar  transfer  functions 


qv  .qv  +  i  q  n  _  7  I  1"^7  Lv 
o  l22  l21U  V-ir  ZoL21  * 

q  -  1,2, 

...  ,a; 

CM 

It 

> 

,V 

2.1 

qj  -L2i<’  -  z„Ln>''‘J  • 

q  -  1,2, 

•  •  • 

j  =  1,2,  ... 

,k 

2.2 

k#  •rk[1*l11(l-20L„)-1Z0]ll[2i 

7T 

II 

r\> 

...  tk^ 

< 

II 

ro 

,V 

2.3 

and 


■  rkL,l(l  -  20Ln)-'cj  ;  k,J  .  1,2,  ... ,k  2.4 


Here,  £  and  denote  the  number  of  external  system  inputs  and  outputs 

ki  2 

which  are  typically  few  in  number.  As  such,  the  e  J  array  composed  of  k^ 
entries  dominates  the  data  base.  Also  note  that  all  of  the  entries  in  the 
data  base  are  formulated  in  terms  of  the  nominal  component  values  and,  as 
such,  the  data  base  may  be  generated  off-line  without  a  priori  knowledge 
of  the  perturbations  to  be  analyzed.  Finally,  the  entire  data  base  may  be 
generated  with  the  aid  of  only  a  single  ji  by  £  (sparse)  matrix  inverse. 

Now,  if  we  assume  that  only  a  single  parameter  is  perturbed,  i.e. 

Z1  =  ck5krk  2.5 

for  some  fixed  k  =  1,2,  ...  ,j<,  to  retrieve  sqv  from  the  data  base  we  must 
evaluate 


:qv  =  | qv 


[ZQ  +  ck5krk]L11)'1[Z0  +  ck5rk]L^2 


k 

in  terms  of  the  elements  of  our  data  base  and  the  variable  parameter,  6  . 
To  this  end  we  invoke  Householder's  formula^ 
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2.7 


(w  +  xy)'1  =  if1  -  w“ 1  x ( l  +  Yw^xrVw'1 
with  VI  =  (1  -  ZqL^),  X  =  -ck6k,  and  Y  =  rkL^  obtaining 
(1  -  [ZQ  +  ck5krk]tlir1  =  C(1  -  Z0Ln)  +  (-ck5k)(rkL11)]'1 


■  I'-Z.L,,)-’  ♦  (1-ZoL1l)'1cks'C(1-4'CLlltl-Z0Ln)'1ckik)'1rkL11(l-Z0L11)-1 


1  k„k  kt 


,-1  .  d-Vul  cVr\„(l-Z0L 

1  1  /  "  “  Is  L/C 


1  -  6kekk 


Now,  upon  substitution  of  2.8  into  2.6  we  obtain 


Sqv  -  L 22  ^  ^in-[Zo-Hck5krk]L11)'1CZ0-Ck6krklL^2 


=  1%  +  1-5,  (1-Z  L,.  )_1[Z  +ck6krk]L^ 


'22  21  v  o  il '  LV 

d  n  -t  i  \  - 1  _k  ,k  k,  / 


•In  j__k.k_kni  v 


L2i(1-zoLn^  c  5  r  Lii(1‘z-L”)'  [z-+c  6  r  ]L 


•nvl  yn'  Lv 

i  kk 

1  -  5  e 


*Mkk,  ,,  ,  ,  v-1,  ,  v.  /  ,k>2Kqk  kk  k,  v 

av  ,  k.akk.v  .  «  *>  r  l11<1-ZoLT1)  ZoL12  +  (i  1  b  e  r  L1 

=  so  *  5  b  r  L12  *  - -  .  {kskk 


■*r 


.k.ak.kv  / ,  k  \  2  r  Kqk  kk  k,  v  Kqk_  kck  k,  v  ,  .  .  . 

!  b  -  d  +  (5  )  [-tr  e  r  L12  +  b1  e  r  L12j  qv  .  5kbqkdkv 

+  - - — r-; - -  —  =  S_  +  r  r 


1  -  5kekk 


,  .k  kk 
1  -  o  e 


which  is  the  desired  symbolic  transfer  function. 

If  we  assume  that  two  parameters  are  oerturbed;  that  is 

2,  -  ck5krk  *  cW  2.10 

a  similar  formula  can  be  obtained  wherein  Householder's  formula  is  applied 
twice.  Since  this  formula  is  subsumed  by  the  general  retreival  formula  derived 
in  the  following  section,  we  simply  state  the  result  without  proof.  In 
particular, 
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III.  Retreival  Theorem 


As  is  apparent  from  equation  2.11,  our  retreival  formulas  are  quite 
complex,  even  for  the  case  p  =  2  and,  as  such,  a  more  compact  notation  is 
required  if  they  are  to  be  tractable.  To  this  end,  we  assume  that 

i< 

i  ;  k  5  1,2,  ...  ,p:  denote  the  potentially  variable  parameters  and  that 

R  \c  V  k 

Z,  =  l  cK6Kr  =  CAR  3.1 

1  k«l 


Of  course,  the  same  expression  applies  to  any  set  of  p  potentially  variable 
parameters  given  an  appropriate  change  of  the  index  set.  To  obtain  the  re¬ 
quired  symbolic  transfer  function  for 


S  .  L22*L2,0-CZo*213Ll,)-1CZ0+Z13L12  3.2 

with  the  above  specified  Z-j  we  now  define  the  following  matrices  made  ud  of 
elements  from  our  data  base 


,11 

12 

,lv 

so 

so  ••• 

so 

s21 

22 

s2v 

so 

so  ••• 

sj'  s f  ...  sj!- 
0  0  0 


0 

5o  •" 

50 

b'1 

b12  ... 

blp 

CM 

-Q 

b22  ... 

b2p 

b^1 

t£2  ... 

b3p 

3.3 


3.4 
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while  A  is  defined  as  per  equation  1.12. 

THEOREM:  Using  the  above  notation 

S  .  42*L210-:Z0Z1]l-n)',CVZ,3L12  -  S0  *  BO  -  iEr’iD 


3.5 


Proof:  First,  we  observe  that 

So  =  L22  +  L21^’ZoL11)  lzoL12 
is  just  the  nominal  system  transfer  function  matrix  while 

8  =  L21(1'ZoLll)  1c 
and 

D  =  RCl  +  l‘11(l-Z0L11)"1Z0]U12  •  R(1-L1,Z0)'1 
via  Householder's  formula.  Finally, 

e  =  Rl]1<i-z0l),)-1c 

where  R  and  C  are  as  defined  by  equations  1.10  and  1.11.  As  such, 
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(1  -  aE)'1  =  (1  -  arl11(i-zol11)‘1c)'1 

»  n  +  ARLn(l  -  (1-ZoL11)"1CaRL11)"1(1-ZoL11)'1C] 

=  [1  +  ARLn(1  -  ZQLn  -  Z^^C] 

=  [1  +  ARL11(1-ZL11)"1C]  3.10 

where  we  have  invoked  Householder's  formula  with  Z  =  1,  X  =  ARL^ ,  and 
Y  =  (1-ZgL.ji  )~^C;  and  equation  1.9.  As  such, 

So  +  B(l  -  AE)_1A0  =  S0  +  L21(1-Z0L11)‘1CC1 +ARL11(l-ZL11)'1C]AR(l-L11Zo)‘1L12 

"  +  ziLn(1‘ZLn)"1zi](1"Lnzo)"lLi2 

=  s0  +  L2i(1'z0Ln)*1{c(1'ZLn)+ziLn](1'ZLn)"1}Zi(1'Liizo)"lLi2 

*  so  +  L2i(1-zoLn)"^1-zoLii)(1-ZLn)"lzi(1-Lnzo)"lLi2 

=  So  +  L21(1"ZL11^  ^l^'HlV  1l12 
’  L22  +  L21(1-ZoLHr1zoL12  +  L21(1'ZLn)'1zi(1'LnZo)"1|-12 

=  L22  +  L21Zo(1"L11Zo)  1|_12  +  L21(1"ZL11^  L12 

=  L22  +  L21[Zq  +  (1-ZL^)  Z1 3(1-L1  -|Z0)  L12 

=  L22  +  L210-ZLn)  C(1-ZL]1  )Zq  +  z-|  ](1-L11Z0)  Li2 

=  L22  +  L21  ( 1  _ZL1 1 )  1[Z  "  ZLnZo]^‘LllV  L1 2 

=  L22  +  L2i(i-ZLii)'lzcl-LnV(1-LnVlLi2 

=  L22  ^  L21 ( 1 "ZL1 1 )  1zL12  =  S 


as  required.  ///////// 
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IV. 


Sensitivity  Formulae 


If  one  is  working  directly  with  the  component  connection  model,  it  is 

4  i 

well  known  that  the  sensitivity  of  S  with  respect  to  a  parameter,  6  ,  can  be 

computed  via  the  formula 

[1  +  L11(l-ZLn)‘1Z]L12  4.1 

and  hence  it  is  appropriate  to  ask  whether  or  not  such  a  sensitivity  matrix 
can  also  be  computed  from  our  data  base.  Since  the  expression 

S  =  S„  +  B(1  -  aE)_1aD  4.2 

o 

is  formally  identical  to  1.4,  if  1  £  i  ^  p  we  may  write 

=  B(1-aE)-1M1[1  +  E(1-aE)_1a]D  4.3 

where 


M. 

i 


dA 

dS1 


L 


0 


0 


1 


0 


4.4 


with  the  one  aDpearing  in  the  ith  diagonal  entry.  Clearly,  the  expression 
can  be  computed  directly  from  the  d-‘ta  base  with  the  same  'evo'  -* 
effort  as  required  for  the  retreiva1  formula. 

In  the  case  where  51  is  no*  ,nc' 
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MAY  -81  R  -  LIU.  K  NAKAJIMA.  P  OLIVIER 


UNCLASSIF IEO 


deviate  from  nominal,  i  >  pin  our  notation,  we  must  first  augment  the  B,  E,  D, 
and  a  matrices  to  include  61  and  then  apply  equation  4.3  to  the  augmented 
system.  To  this  end  we  let 


and 


B1  = 


E1  = 


A3  = 


"b1T 

b12 

...  blp 

b]l'_ 

b21 

b22 

b2p 

b2i 

b*1 

b*2 

b3P 

b & 

d11 

d12 

...  d1^" 

d21 

d22 

...  d2^ 

dp1 

dp2 

...  dp* 

d11 

d12 

...  d1*. 

r  n 

e 

12 

e 

...  elp 

eli_ 

e21 

e22 

...  e2p 

e2i 

epl 

ep2 

epp 

epi 

il 

_e 

i2 

e 

...  eip 

e1  i  _ 

A 

0 

0 

0 

4.5 


4.6 


4.7 


4.8 


The  we  obtain  the  retreival  formulae 

S  =  S  +  B1 (1  -  AaE1)'1AaD1 


and 


[3S‘ 

Idi1 


*  B1(l-AaEi)"1M  [ 1 +E 1 ( 1 -AaE1 )" 1 Aa]D1 

p+1 
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4.9 


4.10 


V.  Updating  the  Data  Base 

In  many  applications  one  uses  a  data  base  such  as  that  described  above, 
as  a  design  tool  to  aid  in  simulating  the  effects  of  various  proposed 
modifications  to  the  system.  When  such  a  modification  is  finally  implemented 
it  is  then  necessary  to  update  the  data  base  to  reflect  the  new  nominal 


parameter  values 


-  P  I,  If  u 

Z„  =  Z  +  T  cK6  r  =  Zn  +  CaR 
0  0  k=l  0 


with  the  aid  of  Householder's  formula  we  may  compute 

n-z0Ln> -  cn-Zohi*  -  CaRLh]"1  =  n-V-n)*1 

+  (1-Z0L11)-1C[1  -  arl11(i-zol11)-1c]'1arl11(i-zol11)‘1 

»  (1-Z^)*1  +  (l-Z0L11)'1C(l-AE)’1ARLlin-Z0L11J’1  5.2 

which  upon  substitution  into  equation  2.4  yields 


ekj  =  ekj  +  [ekl  ek2  ...  ekp](l-4E)'1 a  [elj 

fl2j 


Simi larly. 


bqj  =  bqk  +  [bql  bq2  ...  bqp ]( 1 -aE ) ~ 1  a 


1,  e 


=  dkv  +  [ekl  ek2 


.  ekp](l-AE)-1A  [dlvl 
j2v  I 
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and 


* 


-  s"v„  ♦  [bql  bq2 

o  o 


bqp](1-AE)*1A 


5.6 


As  such,  the  entries  in  our  data  base  can  be  updated  with  a  computational 
effort  which  is  commensurate  with  that  required  by  the  retreival  formula. 


no 


VI .  Examples 

Consider  the  simple  RC  op-amp  circuit  shown  in  figure  1.  The  component 
connection 


Figure  1:  RC  Op-amp  circuit, 
model  for  this  circuit  takes  the  form 


1C 

0 

<T 

vc 

Vr 

= 

0 

R 

0 

V 

v2_ 

0 

0 

K 

Vl_ 

— 

_ 

0  -1  -1 

*c 

1 

vi 

1  0  0 

vr 

+ 

0 

ij 

0  1  0 

v2 

1 

3  ‘ 

•  ] 

^c 

+ 

?] 

Thus  if  all  components  taken  to  have  nominal  values  of  1  we  obtain 


6.1 


6.2 


6.3 


1  s 

-1  1 

0  1 


s 

0 

1 
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6.4 


o-z 


1  0  -S 

1  1  -s 

-1  -1  s+1 


O-z0Lnr  zo 


hi(1-zoLnr 


s  0  -s 

s  1  -S 

s  -1  s+1 


0  0-1 
1  0  -s 

•1  -1  $ 


Ln  +  Ln(1‘zoLi^  zo 


i  o  -i 

s  1  -s 

s  -1  S+l 


,,  we  represent  perturbations  In  the  parameters  C,  9,  and  K  via  the 


matrices 


si  «'  0  0  5 


c'sV  ■  o 


ol  a2  [o  '  5 


cVr2  ■  1 


t 


t 


I 


t 


cVr3 


6.11 


.  k  •? 

Combining  the  appropriate  c  and  r  matrices  with  the  above  expressions  as 
per  equations  2.1  thru  2.4  we  obtain  the  data  base 


s0  *  1  6-12 


b1  =-s 

b2.-i 

b3  *  s+1 

6.13 

d1  -1 

d2  *  0 

d3  =  1 

6.14 

and 

e11  =0 

e,2-0 

e13  *  -1 

e21  *  s 

22 

g*-*-  s  Q 

e23  -  -s 

e  *  -s 

32 

e  =  - 1 

33  c 
e  *  s 

6.15 

where  we  have  deleted  the  q  and  v 

indices 

since  we  are  dealing  with 

a  single- 

input  single-output  system. 

Now,  if  one  desires  to  compute  the  symbolic  transfer  function  with  respect 
to  perturbations  in  the  op-amp  gain  we  have 


*#-  .3,  _  b353d3 

S(s’°  }  so  +  T  7T33 

1-5  e 


1  +  63 
1  -  53S 


6.16 


3 

Recalling  that  5  represents  a  perturbation  from  a  nominal  parameter  value  of 

3  3 

Kq  *  1  our  actual  gain  is  K  a  KQ  +5  *  1+5  ,  which  upon  substitution  into  6.16 

yields 

S<s-K>  ■  i  6-'7 


I 
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which  is  the  classical  gain  formula  for  such  a  circuit. 

Finally,  if  we  desire  to  update  our  data  base  to  reflect  a  new  nominal 
value  for  the  circuit  parameters  of  C  *  1,  R  *  1,  and  K  *  2  we  invoke 

3 

equations  5.3  thru  5.6  with  5*1  yielding 


6.18 

6.19 


and  similarly  for  the  other  elements  of  the  data  base. 


114 


VII .  Conclusions 

The  preceeding  development  has  been  motivated  by  operational  and  main¬ 
tenance  considerations  rather  than  the  design  considerations.  In  such  an 
environment  one  typically  deals  with  a  fixed  nominal  system,  but  carries 
out  repeated  analyses  thereon.  As  such,  the  cost  of  generating  the  required 
data  base  is  secondary  compared  to  the  cost  of  storing  the  data  base  and 
retreiving  information  therefrom.  In  these  respects  we  believe  that  our 
data  base  is  near  optimal.  Since  the  number  of  system  inputs  and  outputs 

p 

is  typically  small  our  data  base  contains  approximately  J<  elements 
2 

(actually  +  kjv+cj.)  +  vgj  where  k.  is  the  total  number  of  parameters 
which  are  potentially  variable.  This  data  base,  however,  contains  sufficient 
information  to  permit  one  to  retreive  symbolic  transfer  functions  for  any 
number  p  ±  k_  of  variable  parameters.  Indeed,  the  number  of  variable  para¬ 
meters  in  a  symbolic  transfer  function  is  reflected  only  in  the  cost  of 
retreival  which  is  on  the  order  of  p  multiplications  (actually  p  +  p  v+ 
pv(g+l)).  Since  p  is  typically  small,  say  five  or  less,  this  is  minimal. 
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Abstract 

This  paper  generalizes  the  classical  Householder's  formula  to  certain 
nonlinear  operators.  This  class  of  nonlinear  operators  is  shown  to  be 
common  in  circuit  theory.  Several  examples  are  provided  that  show  where 
these  operators  occur  and  the  result  is  applied. 
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I.  INTRODUCTION 


The  purpose  of  this  paper  is  to  present  a  technique  for  analyzing 
lumped  analog  systems  with  some  linear  and  some  nonlinear  elements.  It 
is  shown  that  such  a  system  is  described  by  an  operator  of  the  form 

B+YoD  (1.1) 

where  B  and  D  are  linear  and  Y  is  nonlinear.  Since  no  assumptions  are 
made  about  the  nature  of  the  nonlinearities,  it  is  impossible  to  view  the 
operator  YoD  as  small  in  any  sense,  hence,  YoD  has  to  be  viewed  as  a  large 
nonlinear  perturbation  of  the  linear  operator  B. 

The  technique  to  be  presented  is  based  on  a  theorem  that  allows  us  to 
invert  (1.1)  in  two  steps.  First,  invert  the  linear  operator  B.  If  there 
are  NL  linear  elements  and  nonlinear  elements,  B  will  be  an  (NL+NN)x 
(NL+NN)  matrix.  Second,  invert  a  nonlinear  operator  of  rank  N^.  That  such 
a  result  exists,  is  not  surprising.  Those  experienced  in  solving  equations 
involving  such  operators  apply  Gaussian  elimination  until  there  are 
nonlinear  equations  in  unknowns.  Another  way  to  see  that  this  segregation 
can  be  accomplished  is  to  view  the  nonlinear  elements  as  a  "load"  on  an 
appropriate  linear  circuit,  in  much  the  same  way  as  a  circuit  where  one 
nonlinear  element  is  analyzed  by  viewing  that  element  as  the  load  and 
finding  the  Thevenin's  Equivalent  circuit  that  it  sees. 

The  main  result  of  this  paper  is  obtained  by  generalizing  to 
operators  the  form  B+YoD,  a  classical  theorem  concerning  linear  operators 
known  as  Householder's  Formula.  This  classical  result  and  its  generalization 
are  stated  and  proven  in  section  2.  In  section  3,  we  show  such  operators  do, 
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indeed,  occur  in  circuit  theory  and  then  two  examples  are  presented.  The 
results  are  summarized  in  section  4. 

Section  2. 

The  classical  Householder's  formula  [1]  provides  a  means  of  calculating 
the  inverse  of  the  matrix  B+CD  in  terms  of  B-1  and  (I+DB-1C)-1.  If  B-1  is 
known  and  if  the  dimensions  of  C  and  Dare  appropriate,  then  a  great  savings 
in  time  and  effort  can  be  realized  using  this  technique. 

Theorem  1 :  (Classical  Householder's  Formula) 

If  B  is  an  NxN  matrix,  C  is  an  NxP  matrix,  and  D  is  a  PXN 
matrix,  then  (B+CD)-1  =  B"1-B"1C(I+DB‘*1C)“1DB"1 . 

In  the  nonlinear  extension,  the  linear  operator  C  is  replaced  by  the 
nonlinear  operator  Y.  The  proof  of  this  extension  looks,  at  first  glance, 
like  the  proof  of  a  linear  rather  than  a  nonlinear  theorem.  To  see  that 
this  is  indeed  a  nonlinear  result,  the  differences  between  the  nonlinear 
and  linear  operator  algebra  will  be  reviewed  by  giving  two  basic  definitions. 

Definition  1:  (Operator  Addition)  Let  f  and  g  be  two  operators  (linear  or 
nonlinear)  with  the  same  domain,  then  the  operator  f+g  is 
by  the  following 

(f+g)(x)=f(x)+g(x)  (2.1) 
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Definition  2:  (Linearity)  an  operator  f  is  linear  if  for  all  x  and  y  in 
its  domain  and  all  scalars  a  and  8 
f(ax+ey)  =  f(ax)+f(ay) 

=af(x)+Bf(y). 

The  argument  distributes  to  the  left  for  all  operators,  but  the  operator 
distributes  to  the  right  only  for  linear  operators.  With  this  distinction 
in  mind,  we  are  ready  to  state  and  prove  our  main  result  which  is  a  closed 
form  expression  for  (B+YD)"^  in  terms  of  B"1  and  (I+DB"^Y)-1  (rf  course, 
operator  multiplication  is  to  be  interpreted  as  composition,  i.e., 

YD=YoD. 

Theorem  2:  If  i)B  and  D  are  linear  operators,  ii)B-1  exists,  iii)y  is  an 
arbitrary  operator  and  iv)B+YD  is  defined,  then 

(B+YD)"1  =  B"1-B"1Y(I+DBY)"1DB“1.  *  (2.3) 

Proof:  Consider  the  operator  X+XYX  where  X  is  linear  and  Y  is  possibly 
nonlinear. 

X+XYX  =  X( I+YX)  =  I+XY)X. 

If  (I+YX)  and  (I+YX)  are  both  invertible  (it  can  be  shown  [2]  that 
one  is  invertible  if  and  only  if  the  other  is)  we  have 

( I+X Y ) - 1 X  =  X( I+YX)-1  (2.4) 

Now  consider  the  identity 

I  =  ( I+YX) ( I+XY)-1  =  I ( I+YX ) - 1 +YX ( I+YX)-1 
=  ( I+YX ) ~ 1 +Y ( I +X Y ) - 1 X 
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where  we  have  used  (4).  Solving  for  (I+YX)'1  yields 
(5)  (I+YX)'1  =  I-Y( I+XY)”1 X . 

Finally,  consider  the  operator  B+YD. 

(B+YD)"1  =  [(I+YDB”1 )B]-1  =  B"1 ( I+YDB'1 )_1 
Letting  X  =  DB'1  in  (5)  yields 

(B+YD)'1  =  B~ 1 C I  ~Y  ( I+DB~  "*  Y )  ~ 1 DB—1  ] 

=  B'1-B'1Y(I+DB”1Y)"1DB~1 . 

To  see  how  this  result  is  useful,  consider  the  case  where  B+YD  is  an 
Nth  order  nonlinear  operator,  D  a  linear  operator  that  maps  | RN-*-| RP,P<N 
and  Y  a  nonlinear  operator  that  maps  |R  -*-|R  .  This  result  allows  the 
solution  of  N  nonlinear  equation  in  N  unknowns  to  be  replaced  by  the 
solution  of  N  linear  equations  in  N  unknowns  and  also  the  solution  of  P 
(recall  P<N)  nonlinear  equations  in  P  unknowns.  Thus,  we  have,  via  a 
closed  form  expression,  ordered  our  equations  and  unknowns  properly  to 
make  maximum  use  of  linear  techniques  and  minimum  use  of  nonlinear  techniques. 

It  should  be  noted  that  the  proof  of  Theorem  2  relied  on  the  fact  that 
B  and  D  were  linear  operators  and  allowed  Y  to  be  arbitrary.  B  and  D  were 
not  assumed  to  be  matrices  and  Y  was  not  assumed  to  map  |RP-*-|RN.  Any  or 
all  of  the  operators  could  be  differential  operators  and  the  result  would 
still  be  valid.  Regardless  of  whether  the  operators  aredifferential  or 
functional,  we  have  succeeded  in  breaking  it  up  into  a  linear  portion  and 
a  nonlinear  portion.  If  there  are  few  nonlinear  components  in  comparison 
with  the  number  of  linear  components,  the  nonlinearities  can  be  viewed  as 
a  perturbation  on  the  linear  system. 

Section  3 

The  purpose  of  this  section  is  to  show  that  operators  of  the  form 
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B+YD  occur  in  circuit  analysis  problems  and  to  apply  Theorem  2  to  two 
examples. 

This  type  of  operator  arises  naturally  in  nonlinear  network  analysis. 
Consider  the  Node  analysis  of  a  network  with  reduced  incidence  matrix  A, 
T3].  Kirchoff's  Laws  are 

(KCL)  Ai  =  0 
(KVL)  v  =  ATe. 

The  branch  equations  might  be 


Where 


1  *  Gv+j^-Gv^+f (v) 

j-  the  branch  current  vector; 
v_-  the  branch  voltage  vector; 

the  current  source  vector; 
v^-  the  voltage  source  vector; 

the  mode-to-datum  voltage  vector; 


G  is  assumed  to  be  an  "invertible"  matrix  of  differential  operators,  f  is  a 
nonlinear  differential  operator,  and  all  branches  are  voltage  controlled. 


If 


i^AGVs-Ais 


then  Kirchoff's  Laws  and  the  branch  equations  can  be  combined  to  yield 


(AGAT)e+Af(ATe)  =  1_. 


(3.1) 


Letting  AGAT=B,  and  Af(-)=Y,  we  see  that  this  operator  is  of  the  desired 
form.  The  typical  situation  is  for  f  to  be  a  function  of  only  a  few  (p<n) 
linear  combinations  of  the  components  of  v^  then  (3.1)  can  be  rewritten  in 


the  form 
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B+Y(CAT)  e=i 


(3.2) 


which  is  precisely  the  type  of  operator  that  is  amenable  to  the  results  of 
Theorem  2. 

We  now  apply  Theorem  2  to  the  problem  of  solving  two  nonlinear 
simultaneous  equations  in  two  unknowns.  Consider  the  following  nonlinear 
equations 


f(X)  =  z 


-ijX,-  1X,+X 


3 

1 


In  order  to  apply  Theorem  2,  f(x)  must  be  put  into  the  form 


(B+YD) (X) 


where  B  is  an  invertible  matrix.  One  way  to  do  this  is 

H  -TT  X, 

f(X)  =  (B+YD)x  =  »5 


1  1 


V1 


+Y([1,0]X) 


where 


Y( . )  = 


0 


Theorem  2  says  that 

X  =  B‘1z-B"1Y(I+OB"1Y)'1DB"1z  =  XL-XNL 

where 


1 

1" 

2*1 

~Z7 

J 

1_ 

3. 

-21 

XNL  *  B'1Y(1+DB'1Y)‘1DB'1Z 
=  B'1Y(1+DB"1Y)_1DBL 
=  B"1 Y( I+DB~^Y)"^27. 
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Now 


{ I+DB"1 Y)"1 27  =  u 


is  equivalent  to 

( I+DB_1Y)y  =  27. 
(I+DB_1Y)u  =  u+Cl.O] 


The  reason  for  choosing  a  functional  example  is  that  for  large  circuits 
or  systems,  the  differential  equations  are  solved  numerically  so  at  each 
iteration  an  operator  of  the  form  B+YD  must  be  inverted.  To  see  that  this 
is  indeed  the  case,  consider  discretizing  the  differential  equations  obtained 
from  the  component  connection  model  [2]  of  a  system.  The  component  equations 
are  assumed  (here)  to  be  given  in  state  form 

X  =  f(X,a) 
b  =  g(X,a). 

where  a  is  the  vector  of  component  inputs,  b  is  the  vector  of  component  out¬ 
puts,  and  x  is  the  state  vector  of  the  components.  The  connection  equations 
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1 


(KVL  and  KCL  equations)  are  given  by 


a  =  L11  L12 

y  L21  L22 


where  u  is  the  vector  of  system  inputs  and  y  is  the  vector  of  system  outputs. 
If  we  order  the  entries  in  all  of  the  vectors  correctly,  we  can  partition 
the  vectors  in  the  following  manner 


aNl  N~ 

"  aLJ  *  b  =  bLJ  *  a"d  *  =  JL_ 


where  the  superscript  N(L)  denotes  entries  associated  with  the  nonlinear 
(linear)  components.  The  discretized  equations  that  the  computer  is  to 
solve  have  the  form 

J  ,  yN  -N,yN  N, 

,x0  dtxk-i  f  (xk’ak>- 

j0  Vk-rflxk+Bai;- 

bkv<xkN'akK>- 

bkL-cxkLtDak- 


a  N,  NN.  N,.  NL.  L+I  N 
ak  _L1 1 bk  L1 1 bk  L1 2uk  * 


L  .  LN.  N  .  LL.  L+|  L  „ 
ak~Lllbk  Lllbk  L12uk’ 

YK=L21 bk+L22bk+L22uk  * 


The  last  equation  is  just  the  output  equation  and  is  not  used  during  the 
iterations.  These  equations  can  be  put  in  the  following  form 
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TWk+Y(DWk)  = 


-Ido 

0 

0 

0 

-1 

0 

1 - 

X 

0 

0 

0 

0 

0 

0 

aN 

ak 

0 

0 

A-dol 

B 

0 

0 

A 

0 

0 

C 

D 

0 

-I 

A 

0 

I 

0 

0 

.  NN 

~Ln 

-lnl 

L11 

kn 

bk 

0 

0 

0 

I 

.  NN 

~Ln 

-lll 

L11 

r,N,vN  aNxn 

1  <v*k> 

rrVkM 

i*i  1  r  1 

9M0C«.a”) 

0 

r  0 

,.,d^ 

0 

0 

0 

0 

f 

O 

1 _ 

— , 

o 

_ 1 

where  D=  1,0  ,  and  I  is  conformable  with 


4.  Conclusion 

The  classical  Householder's  Formula  has  been  generalized  to  certain 
nonlinear  operators.  It  was  shown  that  these  nonlinear  operators  occur  in 
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circuit  theory,  both  in  the  differential  equations  that  describe  the  circuit 
and  in  the  discretized  equations  that  are  used  in  the  computer  aided 
analysis  of  these  circuits.  It  is  hoped  that  this  result  will  be  as  useful 
a  tool  in  the  fault  analysis  of  nonlinear  circuits  as  the  classical  result 
turned  out  to  be  in  the  fault  analysis  of  linear  circuits. 
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Abstract 

A  fault  analysis  algorithm  appropriate  for  time  varying  and  nonlinear 
systems,  is  developed.  The  algorithm  essentially  constructs  an  observer 
for  a  nonlinear  system  that  is  intimately  related  to  the  system  under  test. 
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INTRODUCTION 

Given  enough  time  and  computing  capability  brute  force  searches  will 

identify  possible  fault  sets.  The  real  problem  in  fault  analysis  is  to 

construct  algorithms  that,  in  some  sense,  locate  the  fault  sets  "efficiently". 

"Efficiently",  in  this  context,  means  that  the  fault  isolation  must  be 

done  relatively  quickly  and  with  limited  on  site  computing.  Such  techniques 

2  3  4 

have  been  developed  to  handle  linear  time  invariant  and  digital  systems.  ’ 

These,  however,  make  heavy  use  of  the  defining  properties  for  these  systems, 
and  do  not  generalize.  The  purpose  of  this  paper  is  to  show  that  an  observer 
for  an  appropriate  nonlinear  differential  equation  can  be  utilized,  on  line, 
to  determine  the  values  of  the  system  parameters.  A  technique,  based  on 
optimal  control  theory,  for  constructing  such  observers  is  also  presented. 

OBSERVERS  AND  FAULT  ANALYSIS 

Consider  testing  a  system  that  is  described  by"  the  nonlinear  state 
equations 

x^j  -  f(x-j,a,u,t) 

y  -  g(x1) 

where  x1  is  the  dynamical  state  vector,  a  is  the  vector  of  parameters  to  be 
estimated  (they  are  assumed  constant  over  the  test  time),  and  u  is  the  input  1 

used  in  the  test  procedure.  If  we  want  to  estimate  a  we  need  to  include  it 
in  the  state  vector,  i.e.,  we  want  to  build  an  observer  for  the  augmented 
differential  equation 

! 

i 

•  », 


I 
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f(x-|  ,a,u,t) 


y  *  g(x]). 

If  it  is  possible  to  build  an  observer  that  will  observe  the  subvector 
a  we  have  solved  the  fault  analysis  problem.  It  would  then  be  necessary  to 
justify  our  solution  in  terms  of  time  and  computation  requirements. 


AN  OBSERVER  DESIGN 


We  chose  to  design  an  observer  with  the  following  structure 

A 

+  H(y  -  y)  (H  time  invariant) 


f  (x,a,u,t) 

- 1 

o 

_ 1 

y  =  g(x1). 

we  term  such  an  observer  as  a  Model  reference  linear  time  invariant  observer. 
The  term  "linear  time  invariant"  is  used  because  the  residuals  enter  in  a 
linear  time  invariant  fashion.  The  problem  is  now  that  of  choosing  H.  To 
avoid  involved  stability  considerations  (at  least  initially)  we  choose  H  so 
that  it  minimizes  the  following  function 

^1 

J(H)  =  /  [(x,-x)2+(a-a)23  dt 
*0 

and  hope  that  the  stability  takes  care  of  itself.  The  construction  of  H 
can  be  done  by  solving  the  following  optimal  control  problem^ 
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min  J(H) 


H  F 


2 


XTQXdt , 


-I 

I 


subject  to  the  differential  equations  constraints 


X  = 


x 

a 

x 

a 


f(x1 ,a,u,t) 

0 

0 

+ 

0 

f(x. ,a,u,t) 

H 

0 

(y-y) 


X(tQ)  =  Cx1(t0)Ia(t0)T,x1(t0)T,a(t0)T] 

Note.  1)  H  will  be  dependent  on  the  X(tQ)  used  in  its  construction,  so  when 
it  is  used  to  estimate  a  (when  a  differs  from  a(tQ))  it  has  little  chance 
of  being  the  optimal  H.  So  even  though  we  use  optimization  techniques  to 
construct  H,  it  will  not,  in  general,  be  optimal.  2)  Several  observers  may 
need  to  be  constructed,  each  one  convergent  for  a  in  a  different  region. 

Experience  indicates  that  only  a  few  components  fail  at  a  time.  Because 
of  this  a  reasonable  approach  is  to  construct  an  observer  for  each  component, 
(thereby  minimizing  the  dimension  of  the  augmented  state  vector)  and  estimate 
the  parameters  for  each  component  in  parallel.  Observers  can  also  be  built  for 
the  common  two  and  three  element  faults. 
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CHAPTER  I 


INTRODUCTION 

The  advent  of  microelectronics,  the  ever-increasing  complexity  and 
compactness  of  electronic  circuits,  together  with  a  need  for  higher 
reliability  in  space,  military,  and  even  commercial  projects,  has 
brought  new  problems  to  industry;  test  and  diagnosis  of  electronic 
circuits  is  one  of  them. 

Presently,  a  printed  circuit  board  may  include  several  hundred 
components;  modular  construction  as  well  as  small  geometrical  dimensions 
make  impractical,  even  infeasible,  and  certainly  uneconomical,  ke 
conventional  test  methods  based  on  classical  laboratory  equipment  such 
as  signal  generators,  meters,  oscilloscopes,  and  probes. 

In  either  the  analog  or  digital  case,  fault  detection  and  location 
in  electronic  systems  is  generally  performed  via  measurements  at  a 
limited  number  of  input  and  output  connections.  These  measurements  are 
then  executed  by  computer  test  programs  to  provide  diagnosis. 

Until  now,  algorithms  for  automatically  generating  test  programs 
have  been  concerned  mainly  with  digital  circuits.  Analog  circuits,  on 
the  other  hand,  have  received  far  less  attention,  due  to  several 
reasons:  Analog  systems  are  frequently  non-linear,  and  the  values  of 
the  parameters  of  the  elements  exhibit  large  deviations  [6];  analog 
signals  are  inherently  more  complex  than  digital  signals.  They  occur 
continuously  in  time,  rather  than  at  discrete  times,  and  their  values 
have  infinite  resolution,  instead  of  being  truncated  into  a  finite 
number  of  bits;  most  importantly,  digital  automatic  test  generation  has 
been  successful  due  to  the  simplified  modeling  at  the  logic  gate  or 
higher  level,  rather  than  the  internal  parameter  level  as  in  the  case 
of  analog  systems.  As  a  result,  most  analog  automatic  test  generation 
and  fault  isolation  techniques  require  a  large  computational  capability 
on  the  ATE  or  off-line  computers  [5]. 
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Several  efforts  have  been  made  to  attack  the  fault  diagnosis 
problem  in  analog  circuits.  The  mul ti -frequency  technique  for  fault 
analysis  in  general  linear  dynamical  systems  was  developed  by  N.  Sen 
and  R.  Saeks  [l]  »  [7]  ,  [8]  ,  [9]  ,  was  considered  to  be  more  efficient  and 
advanced  in  terms  of  output  selection  and  reduction.  By  varying  the 
test  frequency  at  the  same  test  points,  the  number  of  test  points  can 
be  reduced  significantly  compared  to  the  case  of  single  frequency 
measurements. 

This  technique  has  been  generalized  to  non-linear  analog  systems 
by  linearization  of  the  non-linear  components.  Unfortunately,  the 
linearization  concept  fails  in  many  cases  [10] .  For  example,  consider¬ 
ing  a  single-loop  circuit  consisting  of  a  power  supply,  a  resistor, 
and  a  tunnel  diode.  At  a  given  bias,  the  non-linear  characteristic  of 
the  tunnel  diode  has  a  particular  slope.  With  the  breakdown  of  the 
bias  resistor,  the  linearization  of  the  non-linear  characteristic  of 
the  tunnel  diode  at  the  new  biasing  point  will  be  different,  and  the 
diode  will  appear  to  be  faulty  if  one  ignores  the  fact  that  the  diode 
is  operating  at  a  different  bias. 

Recognizing  the  linearization  problem  in  multi -frequency  test, 
the  fault-diagnosis  of  non-linear  analog  systems  in  the  D.C.  case  was 
studied  by  V.  svanathan  and  A.  Sangiovanni-Vincentelli  2  .and  N. 

David  and  A.N.  Wilson  [3]  at  the  component  parameter  level.  These  tech¬ 
niques  have  these  disadvantages:  the  required  number  of  test  points  is 
more  than  required  by  the  multiple  frequency  technique;  the  tremendous 
amount  of  computer  time  is  required  to  solve  very  complicated  non-linear 
equations.  Since  this  is  a  D.C.  test,  these  techniques  can  only  be 
applied  to  memoryless  systems  {  without  reactive  components  such  as 
capacitors  and  inductors  ). 

The  dual -mode  technique  for  fault-diagnosis  for  non-linear  analog 
systems  is  introduced  as  a  compromise  between  the  above  approaches. 

Mode  1  (A.C.  Test):  utilizing  the  multi -frequency  technique  to 
search  for  the  faulty  linear  components,  then  the  faulty  linear  component 
values  are  calculated.  The  non-linear  components  of  the  circuit  under 
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diagnosis  are  replaced  by  small  signal  A.C.  or  linearization  equivalent 
circuits.  Therefore,  this  A.C.  test  is  performed  at  the  internal 
parameter  level  for  both  linear  and  non-linear  components. 

Mode  2  (D.C.  Test):  this  test  is  used  to  diagnose  a  reduced  non¬ 
linear  problem  after  the  linear  components  are  determined  in  Mode  1. 

The  nominal  values  or  the  faulty  values  given  by  Mode  1  are  used  for 
computations  depending  on  whether  the  linear  components  are  in  tolerance 
or  not.  A  non-linear  device  or  element  is  treated  as  a  blackbox  with 
inputs  and  outputs.  The  task  is  to  find  the  D.C.  inputs  and  outputs  or 
the  operating  conditions  of  the  blackbox,  not  the  internal  parameters 
of  the  non-linear  devices.  This  approach  is  compatible  to  previous 
successful  techniques  in  digital  systems  performed  at  the  logic  gate  le¬ 
vel  or  higher.  The  fault-diagnosis  algorithm  is  shown  in  Figure  1. 

The  Mode  2  formulation  and  examples  are  presented  in  this  thesis 
together  with  a  comparison  with  existing  techniques  in  fault-diagnosis 
for  non-linear  analog  systems.  Finally,  an  experiment  on  a  four-transis¬ 
tor  amplifier  circuit  is  conducted  to  verify  the  theorectical  solutions 
and  to  determine  the  sensitivity  of  the  dual -mode  approach. 
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Fig.  1  Fault-Diagnosis  Algorithm  by  the  Dual -Mode  Technique 
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Fig- 1 


Fault-Diagnosis  Algorithm  by  the  Dual -Mode  Technique 
(Continued) 
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CHAPTER  II 


THEORETICAL  DEVELOPMENT 

The  theoretical  formulation  of  the  D.C.  test  for  the  non-linear 
problem  is  based  on  the  component  connection  equations  W  : 

a  =  Mjjb  + 
b  =  M21b  +  M22u 

Where  u  and  y  represent  the  vectors  of  accessible  inputs  and  outputs 
which  are  available  to  the  test  systems,  a  and  b  represent  the  compo¬ 
nent  input  and  output  vectors,  respectively.  The  relationship  between 
a  and  b  is: 

b  =  Z  a 

Although  the  symbol  Z  is  used,  the  components  are  not  assumed  to  be 
represented  by  an  impedance  matrix.  Indeed,  hybrid  models  are  used 
in  most  of  our  examples. 

Before  the  circuit  is  analyzed  in  the  D.C.  test,  the  capacitors 
and  inductors,  which  are  assumed  to  be  linear  components,  are  replaced 
by  open-circuits  and  short-circuits  respectively. 

The  component  vectors  a  and  b  are  partitioned  into: 


aL 

Y 

a  = 

b  = 

-3n. 

bN 

w 

Where  a^  and  b^  are  the  linear  component  input  and  output  vectors,  and 
a^  and  b^  are  the  non-linear  component  input  and  output  vectors. 
Therefore: 


when  a  linearized  model  is  employed. 
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The  component  connection  matrix  is  partitioned  accordingly: 


r  ^ 
nn 

mln 

M11 

M12 

r  1 

bL 

aN 

— 

Mu 

mnn 

"ll 

M  N 
M12 

bN 

y 

M21 

M  N 

M21 

1 

C\J 

C\J 

2: 

u 

Equation  (1)  and  (2)  can  be  solved  simultanously  to  yield: 


Thus  the  inputs  and  outputs  of  the  non-linear  devices  can  be 
computed  by  a  few  simple  matrix  operations.  These  values  are  checked 
against  the  operating  characteristics  of  the  corresponding  functional 
devices  for  fault  isolation. 

The  coefficient  matrices  of  (3)  and  (4)  can  be  pre-computed  if 
the  linear  components  are  not  faulty.  The  matrix  should  be  changed 
to  incorporate  the  faulty  values  of  the  linear  components,  if  any,  to 
avoid  the  computational  errors  caused  by  the  use  of  nominal  values  of 
the  linear  components. 

For  each  pair  of  non-linear  component  input-output  signals,  one 
test  point  is  required.  For  instance,  a  bipolar  transistor  can  be  modeled 
with  two  (2)  input-output  pairs: 


1 - 

CO 

i _ 

VBE 

aN  = 

bN 

_rc. 

_VCE_ 

Therefore,  two  (2)  measurements  must  be  taken.  Non-linear  integrated 
circuits  can  be  modeled  in  the  same  manner.  The  number  of  test  points 
required  in  this  mode  goes  up  linearly  with  the  number  of  non-linear 
devices  in  systems. 


The  matrix  +  ^21  s''n9u^ar  if 

points  are  not  chosen  properly.  The  selection  of  test  poin 
the  above  matrix  non-singular  will  be  discussed  in  detail 
Appendix. 


the  test 
ts  to  make 
in  the 
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CHAPTER  III 


EXAMPLES 


Fig.  2(a)  Single-Transistor  Amplifier 


I 

! 


Fig.  2(b)  Small  Signal  A.C.  Equivalent  Circuit 
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Multi -Frequency  Test  (Mode  1)  [1],  [12] 


The  small-signal  A.C.  equivalent  circuit  of  the  circuit  of 
Fig.  1(a)  is  shown  in  Fig.  1(b).  The  connection  equations  associated 
with  the  A.C.  or  mode  1  test  are  as  follows  : 
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0 

0 

Here  we  initially  allow  V_  ,  IC1  ,  VRi  ,and  IE  to  be  taken 
as  test  outputs.  The  measure  of  testability  min  is  used  to  extract  a 
reduced  set  of  test  outputs  from  these  options.  According  to  table  I 
two  (2)  is  the  minimum  number  of  test  outputs  in  this  example  ,  which 

x 

suffices  to  yield  min  =  0  (perfect  testability)  or  to  provide  locally 
unique  solutions  for  the  fault-diagnosis  equations  . 

The  test  measurements  are  taken  at  the  two  output  VR^  and  IE 
at  twelve  (12  =  n  -  <£min  )  distinct  frequencies  ,  where  n  is  the  dimen¬ 
sion  of  the  parameter  vector  r  .The  faulty  parameters  can  be  identified 
by  using  the  Householder's  Formula  and  the  optimization  algorithm[13]  . 
Only  the  faulty  parameter  values  for  the  1  inear  components  of  the  circuit 
need  to  be  calculated  to  be  included  in  the  calculations  in  mode  2. 
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Table  1  Measure  of  Testability  for  the 
Amplifier  Circuit  of  Figure  1 


OUTPUT 


IC1 

VRA 

IE 


MIN. 

3 

2 

2 

3 


VIC1 

Vvra 

v0»it 

ICl.VR^ 

ICl, IE 
VR'.IE 


0 

1 

0 

2 

1 

0 


V0  ,  IC1  ,  VR ^  ;  o 

V0  ,  ICl  ,  IE  i  0 

vo  *  vra  •  IE  I  0 

ICl  ,  VR^  ,  IE  jo 

V0  ,  ICl  ,  VR^  ,  IE  i  0 
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Assuming  that  the  capacitors  of  the  amplifier  circuit  are  not 
shorted.  They  are  treated  as  open-circuits  in  D.C.  test,  and  removed 
from  the  circuit  before  Mode  2  analysis  begins.  The  test  point  measure¬ 
ments  are  chosen  at  the  same  test  points  used  in  the  multi -frequency 
test.  The  connection  equations  for  this  circuit  in  D.C.  test  are: 


0  -1  0  0  0  0 

1000-10 
0  0  0  0  1  1 

0  0  0  0  0  1 

0  1-10  0  0 

0  0-1-1  0  0 

0  -1  0  0  0  0 

0  0  0  0  1  1 


The  non-linear  component  input  and  output  vectors  a^  and  b^  can  be 
found  via  (3)  and  (4) : 
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VR1  “ 
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IB1 
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VCE1_ 

IC1 

The  non-linear  inputs  and  outputs  are  determined  via  (3)  and  (4): 
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where  R  =  R1  +  R2  +  R3 

*  *  RE.RC 
K  EE+RC 
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The  component-connection  equations  for  Mode  2  analysis  are: 


VR1 

VR2 

IR3 

IR4 

IR5 

IB1 

VBE2 

VCE1 
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1  + 
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Fig.  5  Video  Amplifier  Circuit 


The  same  circuit  is  used  as  an  example  in  [3].  In  [3]  ten  (10) 
test  points  are  required  besides  the  input  terminals.  The  solutions 
or  the  values  of  the  parameters  are  obtained  by  solving  eighteen  (18) 
non-linear  equations.  On  the  other  hand,  the  dual-mode  technique  requires 
eight  (8)  test  points  besides  the  input  terminals  and  the  solutions 
can  be  obtained  by  straight-forward  matrix  manipulations. 


* 
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After  the  capacitors  are  removed,  the  connection  equations  are: 
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CHAPTER  IV 


COMPARISON  OF  FAULT- DIAGNOSIS  TECHNIQUES  FOR  NON-LINEAR  SYSTEMS 


Referring  to  the  single-transistor  amplifier  circuit  of  Fig.  6(a) 
and  its  equivalent  circuit  of  Fig. 6(b)  using  Ebers-Moll  Model  [11]. 


Fig.  6(a)  Fig.  6(b) 


where  R  =  R1//R2 

IEF  =  K1(exp(x1v1)-1) 

ICF  =  K2(exp(A2v2)-l) 

The  node  equations  at  C  and  B  are  : 

(V+/RC)  -  ajlf-p  +  IEF  =  (vx  -  v2)/RC  (5) 

(V+/R1)  +  £F  “  l^p  "  Vj/R  (6) 

Substitute  I£F  and  ICF  into  (5)  and  (6)  to  obtain: 

V+/RC  -  ajK2(exp(  x2v2)-l)  +  K^expUjVjM)  =  (vx  -  v2)/RC  (7) 

V+/R1  -  aNK1(exp(x1v1)-l)  -  K2(exp(x2v2)-1)  =  Vj/R  (8) 


Suppose  Vj  is  chosen  to  be  the  test  point  measurement,  by 
solving  (7)  and  (8)  simultanously,  v^  can  be  expressed  as  : 

y  =  vj  3  f  (u,r1,r2,rc) 

where  u  is  the  input  voltage, V+  in  this  case  . 
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The  Jacobian  matrix  of  y  with  respect  to  ,  R2>  and  Rc>  which 
is  the  basis  of  V.  Visvanathan  and  A.  Sangiovanni-Vincentelli  formu¬ 
lation,  must  be  computed  numerically.  This  example,  on  the  most 
simple  circuit,  has  thus  shown  the  impractical ity  of  the  above  method. 

Method  of  Nasrollah  David  and  A.  N.  Wilson  [3] 

Referring  to  the  single  transistor  amplifier  circuit  in  [3].  This 
circuit  requires  five  (5)  test  points,  excluding  the  input  terminals. 
Also,  the  parameter  values  can  be  obtained  by  solving  eight  (8)  non¬ 
linear  equations. 


Dual -Mode  Fault-Diagnosis 

Mode  1 :  referring  to  the  single-transistor  amplifier  circuit  in 
[1],  the  set  of  test  point  measurements  V'Ra  and  1^.  yields  the  perfect 
testability  (<5min  =  0)  implying  that  the  fault  diagnosis  equations 
have  locally  unique  solutions. 

Mode  2:  referring  to  the  circuit  in  Figure  1  of  Example  1  in  the 
previous  sections,  the  same  test  points  can  be  used  to  measure  VRg 
and  I£.  In  general  case,  the  set  of  test  points  for  the  two  modes  will 
overlap  each  other. 


Summary 

Computationally,  the  dual-mode  fault-diagnosis  technique  uses  only 
straightforward  matrix  manipulations.  This  is  an  advantage  over  sol¬ 
ving  non-linear  equations.  The  required  number  of  test  points  in  the 
dual-mode  technique  is  much  less  than  that  of  Wilson's  technique. 

Only  one  test  point  is  required  by  Sangiovanni-Vincentell i 's  method 
but  the  trade-off  is  to  solve  a  very  complicated  set  of  non-linear 
equations. 


172 


CHAPTER  V 


FAULT-DIAGNOSIS  ALGORITHM  FOR  TRANSISTORS 

The  calculations  in  Mode  2  of  the  dual-mode  fault-diagnosis 
technique  are  performed  to  provide  the  D.C.  operating  points  of  the 
non-linear  devices.  In  the  case  of  bipolar  transistors,  the  opera¬ 
ting  conditions  can  be  determined  by  the  collector  current  I q ,  the 
base  current  Ig,  the  base-emitter  voltage  Vg£,  and  the  collector- 
emitter  voltage  Vq£.  The  fault-diagnosis  algorithm  for  bipolar  tran¬ 
sistors  is  developed  on  the  basis  of  their  known  operating  charac 


Table  2  Transistor  Operating  Modes 


Modes 

!b 

ic 

VBE 

VCE 

Active 

!b 

«n.!b 

^  .6  v 

Vce  >  VBE 

Saturation 

!b 

V'8 

1  .7  v 

VCE  <  VBE 

Cut-off 

0 

0 

<  .45v 

VCE  %VCC 

where  bn  is  the  nominal  current  transfer  ratio  in  active  mode,  e$ 
is  the  saturated  current  transfer  ratio,  and  VgC  is  the  power  supply 
voltage  connected  to  the  collector. 

The  algorithm  for  field-effect  transistors  is  even  simpler 
because  the  input  or  gate  current  Ig  is  always  zero.  Field-effect 
transistors  are  characterized  by  the  gate  voltage  Vg,  the  drain 
voltage  Vg,  and  the  drain  current  Ig. 

The  fault-diagnosis  algorithm  for  bipolar  transistors  shown 
in  Figure  7,  is  used  to  analyze  the  experimental  results  in  the  fol¬ 
lowing  section. 
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Fig.  7  Fault  Diagnosis  Algorithm  for  N-P-N  Bipolar  Transistors 
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7  Fault  Diagnosis  Algorithm  for  N-P-N  Bipolar  Transistors 
(Continued) 


CHAPTER  VI 


EXPERIMENTAL  RESULTS 

The  video  amplifier  circuit  of  Example  4  was  built  and  tested 
at  nominal  operating  conditions  and  at  intentionally  faulty  condi¬ 
tions.  The  measured  values  of  the  non-linear  devices 'operating  con¬ 
ditions  are  compared  to  those  obtained  by  calculation.  The  computation 
in  this  particular  example  is  simple  enough  to  be  carried  out  by  a 
programmable  hand-held  calculator.  The  computational  error  are  inves¬ 
tigated  experimentally  in  the  following  cases. 

Case  I.  Nominal  Operating  Conditions 

The  components  of  the  circuit,  which  consist  of  four  (4)  tran¬ 
sistors  and  nine  (9)  resistors,  were  carefully  analyzed  before  the 
experiment  was  started.  The  four  transistors  2N2222A  were  checked  on 
a  curve  tracer.  Their  betas  or  current  gains  varied  between  one  hun¬ 
dred  twenty  (120)  to  two  hundred  sixty  (260).  Thus  there  is  more  than 
one  hundred  percent  (100%)  variation  among  the  various  transistors. 
Furthermore,  these  transistors  are  highly  sensitive  to  temperature. 

For  example,  a  transistor  that  carries  twice  the  amount  of  collector 
current  than  another  will  generate  more  heat  and  change  its  charac¬ 
teristics.  The  temperature  sensitivity  has  a  major  effect  on  the  accu¬ 
racy  of  the  test  results.  All  the  resistors  were  within  five  percent 
(5%)  tolerance.  The  measured  values  and  the  nominal  values  of  the 
resistors  are  compared  in  Table  3.  The  experimental  test-point  mea 
surements  for  all  cases  of  the  experiment  is  tabulated  in  Table  4. 

The  experiment  was  performed  in  a  temperature-controlled  environ¬ 
ment  at  twenty  degrees  Celcius  (20®  C).  The  measurements  at  the  test 
points  were  taken  and  used  to  compute  the  transistors'  operating  con¬ 
ditions  with  two  (2)  sets  of  the  resistor  or  linear  component  values. 
The  two  (2)  sets  of  computed  values  for  the  transistors'  operating 
conditions  are  tabulated  in  Table  5. 

The  errors  that  are  produced  by  using  the  manufacturer's  claimed 
values  for  the  resistors  are  less  than  ten  percent  (10%)  for  the  col¬ 
lector  currents,  the  base-emitter  and  the  collector-emitter  voltages. 


Table  3  Comparison  of  Measured  and  Nominal  Values  of  Resistors 


Resistor 

Measured  values 

Claimed  values 
(Nominal ) 

Percent  error 

R1 

1.189  K 

msm 

-.90 

R2 

3.298  K 

-.06 

R3 

5.654  X 

MSB 

.96 

R4 

1.173  K 

1.2  K 

-2.25 

R5 

.316  K 

.33  K 

-4.24 

R6 

.325  K 

.33  K 

-1.50 

R7 

1.010  K 

1.00  K 

-1.00 

R8 

1.461  K 

1.50  K 

-2.60 

R9 

L - 

3.280  K 

3.30  K 

-.61 

Table  4  Test  Point  Measurement  for  All  Cases  of  the  Experiment  at 
V+  =  26.2  Volts  and  V-  =  -28.0  Volts 


Toct  Units  Case  I  Case  II  Case  III  Case  IV 


wm 

volt 

— 

16.04 

28.52 

15.87 

15.86 

I 

volt 

27.31 

25.66 

27.32 

27.32 

VR7 

volt 

6.04 

0.00 

9.827 

5.336 

VR1 

volt 

.0264 

1.611 

.0263 

.0263 

VR9 

volt 

8.83 

0.00 

0.00 

0.00 

VR8 

volt 

8.800 

0.00 

14.32 

0.00 

VR4 

volt 

9.492 

0.00 

15.04 

16.44 

IC4 

mA 

2.665 

0.00 

12.82 

.00866 

k 


Table  5  Transistor  Operating  Conditions  of  Case  I 


Parameter 

Computed  operating  points 
based  on 

Measured  operating 
points(ii) 

%  Error 

(i)  vs.(ii) 

Measured  R 

Claimed  R  (i) 

IB1 

22.2 

22.00 

22.04  uA 

.18 

IB2 

56.76 

70.33 

56.33  uA 

24.90 

IB3 

47.11 

-169.3 

41.75  uA 

-505.5 

IB4 

31.07 

14.76 

28.09  uA 

-47.50 

IC1 

4.808 

4.855 

4.805  mA 

10.40 

IC2 

8.113 

7.748 

8.094  mA 

-4.27 

IC3 

5.978 

6.038 

5.993  mA 

.75 

IC4 

2.665 

2.665 

2.665  mA 

0.00 

VBE1 

.6636 

.6636 

.656  v 

1.16 

VBE2 

.6640 

.6640 

.655  v 

1.30 

VBE3 

.6740 

.7449 

.674  v 

10.50 

VBE4 

.6520 

.6571 

.649  v 

1.17 

VCE1 

10.85 

10.85 

10.82  v 

.28 

VCE2 

16.71 

16.71 

16.70  v 

.06 

VCE3 

11.36 

11.36 

11.36  v 

0.00 

VCE4 

17.37 

17.37 

17.36  v 

.06 

178 


The  errors  are,  however,  worse  for  the  base  currents  due  to  the  fact 
that  the  base  currents  are  so  small.  Considering  the  transistor  Q3, 
the  base-emitter  voltage  is  positive;  therefore,  the  base  current 
cannot  be  negative  or  the  transistor  is  not  faulty.  This  is  thus  a 
case  of  calculation  error. 

Case  II.  A  Bad  Linear  Component  in  the  Circuit 

The  resistor  is  increased  from  3.3K  to  8.9K.  A  set  of  measure¬ 
ments  is  taken  from  the  test  points,  and  used  to  compute  the  transis¬ 
tors'  operating  conditions  with  the  faulty  value  of  taken  or  not 
taken  into  account.  The  results  are  tabulated  in  Table  6.  The  signi¬ 
ficance  is  that  the  transistor  Q2  appears  to  be  faulty  if  the  faulty 
value  of  R^  is  not  taken  into  account.  The  base  current,  which  is 
equal  to  the  collector  current  in  magnitude,  suggests  that  the  base  and 
the  collector  of  Q2  are  shorted  together  while  the  emitter  is  open. 

On  the  other  hand,  if  the  faulty  value  of  R?  is  used  for  computation 
of  the  base  and  collector  current  of  Q2,  these  currents  are  almost  zero 
(0),  indicating  that  the  transistor  Q2  is  operating  at  the  origin  of 
the  characteristic  curve. 

Case  III.  A  Bad  Transistor  in  the  Circuit 

Transistor  Q4  is  replaced  by  a  known  faulty  transistor,  whose 
emitter  is  open-circuited  and  base  to  collector  junction  is  almost 
short-circuited.  This  is  a  comnon  type  of  failure  when  the  emitter  to 
base  junction  is  forward-biased  too  much.  The  base-collector  voltage- 
current  characteristics  of  the  above  transistor  is  shown  in  Figure  8. 
The  transistors'  operating  points  are  calculated  and  tabulated  in 
Table  7. 

Referring  to  the  calculated  operating  points  based  on  the  mea¬ 
sured  resistor  values  of  Table  7,  the  value  of  IB4  equals  in  mag¬ 
nitude  to  IC4,  and  VBE4  is  much  greater  than  point  six  volts  (0.6V). 
These  data  indicate  that  the  emitter-base  junction  is  short-circuited, 
therefore,  the  transistor  Q4  is  faulty. 
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Table  6  Transistor  Operating  Conditions  for  Case  II 


Transistor 

Parameter 

Faulty  value  of 
is  used  for  operating 
point  calculations. 

Faulty  value  of 

not  taken  into 
account. 

Q1 

IB1 

1.355 

1.355  rnA 

IC1 

3.183 

3.183  mA 

VBE1 

.729 

.729  v 

VCE1 

.020 

.020  v 

Q2 

IB2 

.001 

5.464  mA 

IC2 

-.001 

-5.464  mA 

VBE2 

-2.32 

-2.32  v 

VCE2 

26.20 

26.20  v 

Q3 

IB3 

0.0 

0.0  mA 

IC3 

0.0 

0.0  mA 

VBE3 

0.0 

-0.0  v 

VCE3 

26.20 

26.20  v 

Q4 

IB4 

0.0 

0.0  mA 

IC4 

0.0 

0.0  mA 

VBE4 

0.0 

0.0  v 

VCE4 

26.20 

26.20  v 

fRSCHttwO  Pag*  BUMWXOT  FUJttD 


Fig.  8 


Collector-Base  I-V  Characteristics  of  the  Faulty  Transistor 
Used  in  Case  III  and  IV  of  the  Experiment 
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Table  7  Calculated  Operating  Points  and  Fault-Diagnosis 
for  Case  III 


Q 

Parameter 

Calculated 
value  based 
on  measured 
resistance. 

Diagnosis 

Calculated 
value  based 
on  claimed 
resistance 

Diagnosis 

Q1 

IB1 

22.12 

Good 

21.92  uA 

Good 

IC1 

4.81 

B  =  217 

4.86  mA 

B  =  222 

VBE1 

.6540 

.654  v 

VCE1 

11.01 

11.01  v 

Q2 

IB2 

2.15 

At  rest 

-47.56uA 

VBE>.5v  , 

IC2 

73.47 

IB  -  0 

-517.4  uA 

IB.IC  must 

VBE2 

-4.71 

IC  0 

-4.71  v 

be  zero  : 

VCE2 

11.16 

11.16  v 

Calculation 

error 

Q3 

IB3 

71.80 

Good 

-280.3  uA 

VBE  >.6v 

IC3 

9.73 

B  =  136 

9.83  mA 

IB  must>  0 

VBE3 

.6970 

.813  v 

:Cal .error 

VCE3 

2.05 

2.05  v 

Q4 

IB4 

-12.82 

Bad 

-12.82  mA 

Bad 

IC4 

12.82 

IB=-IC 

12.82  mA 

IB=  -IC 

VBE4 

19.21 

VBE».6v 

19.27  v 

VBE».6v 

VCE4 

26.20 

:B-E  opened 

26.20  v 

C-B  short. 
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Again,  using  the  nominal  values  of  the  resistors  for  calculation 
produces  significant  errors.  Consider  the  transistor  Qg,  the  base  to 
emitter  voltage  is  negative,  therefore,  the  transistor  is  operating 
in  the  cut-off  region  or  both  the  base  and  the  collector  currents  have 
to  be  close  to  zero  (0).  However,  the  calculations  show  a  relatively 
large  negative  value  for  both  of  these  currents.  The  error  also  occurs 
in  the  base  current  calculation  for  transistor  Q3,  which  causes  an 
ambiguous  state  for  the  above  transistor. 

Case  IV.  Transistors  Q3  and  Q4  are  Faulty 

Transistors  Qj  and  are  replaced  by  known  faulty  transistors, 
whose  base-collector  characteristics  are  shown  in  Figure  8.  The  tran¬ 
sistors'  operating  points  are  calculated  and  tabulated  in  Table  8. 

The  calculation  of  operating  points  based  on  the  measured  resistor 
values  is  accurate,  while  those  based  on  the  nominal  resistor  values 
creates  errors. 
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Table  8  The  Calculated  Operating  Points  and  Fault- 
Diagnosis  for  Case  iy 


Q 

Parameter 

Calculated 
value  based 
on  measured 
resistor 
values 

Diagnosis 

Calculated 
value  based 
on  claimed 
resistor 
val ues 

Diagnosis 

Q1 

B 1 

22.12 

Good 

21.92  uA 

Good 

4.81 

8  =  217 

4.86  mA 

8  =  222 

.654 

.654  v 

■ 

11.02 

11.02  v 

Q2 

IB2 

-1.19 

At  rest 

-50.9  uA 

At  rest 

IC2 

74.22 

IB, 1C  2j0 

-244. 3uA 

IB.IC  2-  0 

VBE2 

-6.10 

VBE  <0 

-6.1  v 

VBE  <0 

VCE2 

9.76 

9.76  v 

Q3 

IB3 

-5.28 

Bad 

-5.34  mA 

Bad 

IC3 

5.28 

IB*  -IC 

5.34  mA 

IB=- IC 

VBE3 

18.11 

VBE».6v 

18.20  v 

VBE».6v 

VCE3 

20.86 

20.86  v 

Q4 

IB4 

-8.66 

Bad 

-8.66  mA 

Bad 

IC4 

8.66 

same  as 

8.66  mA 

same  as 

VBE4 

19.25 

Q3 

19.30  v 

Q3 

VCE4 

26.20 

26.20  v 

i 
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CHAPTER  VII 


CONCLUSION 

The  dual-mode  fault-diagnosis  technique  for  non-linear  systems 
has  been  formulated  in  the  preceding  sections  via  the  component- 
connection  equations.  Mode  2  or  non-linear  D.C.  analysis  is  performed 
at  the  device  or  element  level  rather  than  at  the  internal  parameter 
level  of  the  corresponding  non-linear  devices,  which  is  not  only  com¬ 
patible  to  previous  successful  fault-diagnosis  techniques  in  digital 
systems,  but  also  very  practical  in  today’s  increasingly  complex  elec¬ 
tronic  systems. 

Computation-wise,  this  technique  is  much  more  advantageous  than 
the  other  existing  techniques  in  analog  non-linear  systems  because  of 
the  use  of  linear  matrix  manipulations  rather  than  solving  complex 
non-linear  equations.  However,  this  technique  requires  a  relatively 
larger  number  of  test  points  compared  to  the  method  of  V.  Visvanathan 
and  A.  Sangiovani-Vincentelli  [2]. 

The  experimental  results  indicate  that  the  D.C.  or  non-linear 
fault-diagnosis  techniques  cannot  be  used  to  diagnose  the  tolerance 
or  soft-fault  problems  due  to  the  significant  errors  introduced  by 
the  use  of  the  manufacturer's  claimed  component  values.  However, these 
errors  can  be  eliminated  by  the  dual -mode  technique  since  the  measured 
values  of  the  linear  components  can  be  computed  in  Mode  1  or  the  multi - 
frequency  test.  Furthermore,  for  the  Mode  2  analysis,  the  bad  value  of 
a  faulty  linear  component  has  to  be  calculated  so  that  the  more  accu¬ 
rate  results  and  fault-diagnosis  can  be  obtained. 
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APPENDIX 


TEST  POINT  SELECTIONS 
Analysis  of  the  Singularity  of  the  Matrix 

F  =  (  M2Ll(ZL1-Mli)’lMil  +  M21  > 


This  analysis  is  intended  to  expose  the  readers  to  the  problem 
of  the  test  point  selection  in  the  D.C.  analysis  of  the  dual-mode 
fault-diagnosis  technique.  The  matrix  F  can  be  arranged  in  the 
following  form: 


F  -  (  M2Li(Z[1-M^)-1M^  +  ) 


where: 


I  =  identity  matrix 


The  matrix  depends  totally  on  output  selections;  therefore, 
the  singularity  of  matrix  F  is  also  determined  by  the  test  point  or 
output  selections. 
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Referring  to  Example  2  of  Section  III,  the  matrix  K1  is: 


R2+  R3 


R 

re 

R 

_re 

1 

1 

V  rc 

re  +  Rc 

re  +  rc 

-R1R2 

R2R3 

0 

0 

R 

R 

-R1R3 

R3^R1+  R2^ 

0 

0 

R 

R 

• 

rerc 

rerc 

-RE 

-RE 

re  *  Rc 

re  +  Rc 

re  +  Rc 

RE  +  RC 

9 


If  V^,  IR1,  IRC,  VR3  are  chosen  as  test  measurements,  where 
=  +  VCE 2  the  matrix  is: 


0 

0 

0 

0 

1 

0 

0 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

The  resulting  matrix  F  is  singular  because  row  1  and  3  are 
linear  dependent. 


ReRc 
Re+  Rc 

R2+R3 

R 


RC+ 


_RjR3 


Re+  Rc 


Rc 

re+rc 


Rc 

re+  rc 


Rq+  re 


-R3(rI+  R2) 

R 


rc+  re 


-1 

Rc+  Re 
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However,  with  a  minor  change  in  output  selection,  V^,  IR1, 
IRC,  VR3  are  chosen,  the  first  row  of  becomes: 

C°  0  0  0  1  jo  0  0  3  , 


and  the  first  row  of  matrix  F  becomes: 


rerc 

re+rc 


rerc 


-Rr 


rE+  Rc  Re+  Rc  Re+  Rc 


Therefore,  F  is  non-singular. 

Referring  to  the  Example  of  Section  III,  the  matrix  K*  is: 


1/R1  000 

1/R2  000 


0  R3  R3  0 

0  0  0  R4 

0  R5  0  R5 


1 


If  VM>  I+,  VB1,  VB2  are  chosen  as  test  points,  the  matrices 


M21  and  F  are: 


'21 


■» 

0 

0 

0 

-1 

0 

1 

1 

1 

0 

0 

0 

-■ 

0 

1 

0 

0 

0 

0 

1 

1 

1 

0 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

•» 

0 

-1 

0 

0 

1 

1 

1 

0 

0 

0 

0  _ 

F  = 


0  0  0  -fy 

-1/R1  1  1  1 

10  0  0 

0  -R3  -R3  0 


The  F  matrix  is  singular  because  the  columns  2  and  3  are 
linearly  dependent. 


However,  if  VM>  It,  VR5,  VB2  are  chosen,  the  matrix  F  is 


non-singular: 


"21- 


0 

0 

0 

-1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

1 

1 

0 

0 

0 

0 

F  = 


■1/R1 

■1/R1 

0 

0 


0 

1 

R5 

-R3 


0 

1 

0 

■R3 


0 

1 

R5 

0 
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Considering  the  second  and  third  column  of  matrix  K*,  the 

differences  between  these  two  columns  are  the  elements  kA  =  Rc, 
111  52  b 

K&3  *  0  and  =  1,  Kg3  *  1.  A  subset  of  Vl^,  which  consists  of 

two  outputs,  has  to  be  chosen,  in  such  a  way  that  its  product 

with  the  second  and  third  column  of  K*  is  non-singular.  This 

condition  is  satisfied  when  one  of  the  outputs  contains  VR5, 

as  in  the  third  row  of  of  the  preceding  discussion,  which 

results  a  non-zero  element  at  row  3  and  column  2  of  matrix  F. 

Referring  to  Example  4  of  Section  III,  the  matrix  K*  is: 


In  order  to  make  det 


0  R9  R9 
R4  -R4  0 
1  0  1 


0,  one  of  the  elements  of 


the  third  row  of  this  matrix  has  to  be  changed  to  aero  as  follows 


- 

0 

R9 

R9 

0 

R9 

R9 

R4 

-R4 

0 

or 

R4 

-R4 

0 

0 

— 

0 

1  _ 

1 

0 

0 

.  „  t,e.e-  cttiprt  the  last  output  measure 

Physically,  this  means  one  has  to  select  tne  « 

ment  such  that  it  does  not  contain  both  IC4  and  IC3. 
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